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Abstract 

"^ ' This is the second in a series of two papers. While in Paper I [l| we derive semiclassical Boltzmann 

\^ ' transport equations and study their flow terms, here we make use of the results from that paper 

2 ' and address the collision terms. We use a model Lagrangean, in which fermions couple to scalars 

O through Yukawa interactions and approximate the self-energies by the one-loop expressions. This 

r^ '■ approximation already contains important aspects of thermalization and scatterings required for 

I [ quantitative studies of transport in plasmas. We compute the CP-violating contributions to both 

qj [ the scalar and the fermionic collision term. CP-violating sources appear at order h in gradient and a 

'-^ I weak coupling expansion, as a consequence of a CP-violating split in the fermionic dispersion relation. 

. ^ [ These represent the first controlled calculations of CP-violating collisional sources, thus establishing 

/S ■ the existence of 'spontaneous' baryogenesis sources in the collision term. The calculation is performed 

d by the use of the spin quasiparticle states established in Paper I. Next, we present the relevant leading 

order calculation of the relaxation rates for the spin states and relate them to the more standard 

relaxation rates for the helicity states. We also analyze the associated fluid equations, and make a 

comparison between the source arising from the semiclassical force in the flow term and the collisional 

source, and observe that the semiclassical source tends to be more important in the limit of inefficient 

diffusion/transport . 
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1. INTRODUCTION 

This is the second part in a series of two closely related papers, where we study the emergence of 
Boltzmann transport equations from quantum field theory for particles with a space-time dependent 
mass, which often arises via a coupling to a varying background field. This setting is typically 
realised at phase transitions, an important example being electroweak scale baryogenesis, in which 
case particles acquire mass at the electroweak phase transition through the Higgs mechanism. This 
paper is an immediate continuation of Paper I [l|, to which we refer for a general introduction and 
motivation. We frequently refer to sections and equations of Paper I, and indicate this by adding a 
prefix "I." to the respective reference. 

We work in a model (1.2. 13-1. 2. 14) with chiral fermions ipL,R that couple to scalars via a Yukawa 
interaction Ci^t (see Eq. ()2.H) below): 

C = t^^^- ^luiiPr - ^Rm*^L + (9^0)^ {d^<P) - (t^'^M^ + £i„t . (1.1) 

Due to the interaction with some background field the masses m and M^ are space-time dependent, 
and they may contain complex elements to allow for CP-violation. In Paper I, based on a 2PI 
effective action we derive the Schwinger- Dyson equations for this system, from which we arrive at 
the Kadanoff-Baym equations (1. 2. 54-1. 2. 55) for the scalar and fermionic Wightman functions A^'^ 
and S^'^ , which in the Wigner representation read 

[e ~-d^ + ik-d- MV^^-^^) A<'> - e-**{n4{A<'>} - e-'^{n<'>}{AJ = C^ (1.2) 
( ^ +^ ^ -m,e-5^-«'= - ^75^,e-^^-^'=)s<'> - e-^^{S,}{5<'>} - e'^^{S<'>}{S,} = C^ , (1.3) 

where n*^'^, II/j (S^'^, S^) denote the scalar (fermionic) self-energies. The scalar and fermionic 
collision terms are of the form 

c, = V^^({n>}{A<} - {n<}{A>}) (1.4) 

C^ = ie-^({S>}{5<}-{S<}{5>}), (1.5) 

where o{a}{&} = {l/2){dxa ■ d^b — d^a ■ dj)). In section 1.3 we argue that the last two terms on the 
left hand side of the equations, which represent coUisional broadening due to multiparticle effects, 
can be neglected. 

We then show that to first order in gradients (equivalently, to linear order in K) the scalar Wigner 
function is solved by the classical on-shell distribution function (1.4.12) 

i/\<{k, x) = — [5{k^ - uj^) - 6iko + Lu^)] n^{k, x) , (1.6) 

where the dispersion relation is given by uj^ = (fc^ -|- M^)^/^. This solution applies to the one 
field case, but also to the diagonal entries of the Wightman function in the mass eigenbasis in 
the case with mixing scalars, provided no eigenvalues Mf of M^ are nearly degenerate, that is 
hk ■ d -^ Mf — Mj {\/i,j). From ()1.6|) it then follows that particle and antiparticle densities are the 
on-shell projections (1.4.31, 1.4.34) 



ft{k,x) = n'^{u^ik,x),k,x) (1.7) 

ft{k,x) = n^''P{uj^,k,x)=-[l+n^{-uj^,-k,x)]. (1.8) 

When working to order h accuracy, one gets the following Boltzmann equation for the CP-violating 
density, 6f^ = ft - ft (1.4.36) 

Q^^tX^ (AMMIqA sfik^ ^) = r^ [c^^k, x) + Cl{k, x) + C^{-k, x) + C;(-A:, x)] . (1.9) 
uj^ Zu^ J Jq 2n ^ 

The discussion of the fermionic kinetic equation in Paper I begins with the observation (see Refs. J3| 
and [3|) that, if masses only depend on the z-coordinate {e.g. caused by planar phase interfaces), the 
following spin operator (1.5.18) 

^.(A:) = i(fco7°-^r7)7Y (1-10) 

ko \ / 

is a conserved quantity. Here (cq = sign(/co)(fco — kl — kyY^'^ is obtained by boosting fco to a frame, 
in which k^ and ky vanish. With the spin projector Ps{k) = (1 + sSz{k))/2 {s = ±1) the fermionic 
Wigner function can accordingly be cast to a spin diagonal form (1.5.25) 

Ss = iPs [si'i'g'o - s^'gl + Igt - t^'g'^] ■ (l-H) 

In stationary situations the densities g^ {a = 0,1, 2, 3) are (real) functions of z and momentum. Only 
one of these functions is independent. Indeed, by making use of the constraint equations, one can 
express gl2 3 ^^ terms of gQ. 

In Paper I we show that, when the coUisional broadening (width) is neglected, quasiparticle picture 
remains valid to order h (see also J2] and [2|), such that the constraint equation for the vector density 
g^ is solved by the following on-shell form (1.5.71) 

grik,x) = 2nJ2^-^p^^^\h\n., (1.12) 

where Ug = ns{k, x) is a distribution function, to be determined by solving the kinetic equation, and 
Zg is the rate of change of the energy-squared around the poles (divided by 2uJs). The dispersion 
relation ojg = ujs{k, z) acquires a CP violating, spin dependent, shift due to the interaction with the 
background field (1.5.72): 

uJs=ujq--z \m\'^dze . (1.13) 

Here 9 is the phase of the complex mass and ujq = [k"^ + |mp)2, uq = [k1 + |mp)2. The densities of 
particles and anti-particles with spin s are the on-shell projections (1.5.120-121) 

fs+{k,x) = ns{ujs,k,x) (1.14) 

fs^{k,x) = nf{ujs,k,x) = l-n^s{-i^s,-k,x). (1-15) 



Besides an equilibrium part and a CP-even deviation from the equilibrium density, these particle 
densities contain spin dependent contributions 6fs±- From these we form two combinations, which 
are related to vector and axial density (1.5.131-132): 

6f: ^ 6U-6fs- (1.16) 

Sf: ^ SU - 5f-s- ■ (1.17) 

While the Boltzmann equation for the vector part does not have a semiclassical source (1.5.133) 

r°° dk 
{dt + vo-V, + F,d,:)5f:= ^{^(k)-ICrik)), (1.18) 

Jo ^ 

there is a CP- violating semiclassical source sS^""" (1.5.135) in the equation for the distribution func- 
tion related to the axial density (1.5.134): 

{dt + vo ■ V, + Fodk^)5f: + s5«°- = / ^ (mk) - JCr'ik)) , (1.19) 

Jo '"" 

where Vq = k/utQ and Fq = —{dz\'m\'^)/{2ujQ). The collisional contributions on the right hand side are 
essentially the traces of the spin projected fermionic collision term. 

In Paper I we formally kept the contributions of the collision term in all expressions, as can be seen 
in Eqs. ()1.9p and p.l8ll.l9|) . but we never evaluated them. This will be done here. In the actual 
calculations we truncate the self-energies in the collision term at the one-loop order. 

Section |21 is devoted to a calculation of CP-violating sources both in the scalar and in the fermionic 
collision term. We emphasize that this represents the first ever calculation of CP-violating sources 
arising from nonlocal (loop) effects. It would be incorrect to associate or identify our calculation 
with, for example, Ref. |^, where the author evaluates the self-energy for stops and charginos by 
taking account of the contributions from two mass insertions, which represent an approximation of 
the stop and charginos interactions with the Higgs background. Formally, all mass insertions are to 
be considered as singular contributions to the self-energy, and hence ought to be viewed as part of the 
(renormalized) mass term. We emphasize that in the formalism presented here all contributions to 
the self-energy arising from the mass insertions (resummed) are fully included in the fermion masses 
in the flow term. Treating mass insertions in any other way is not in the spirit of Kadanoff-Baym 
equations, and hence it is not guaranteed to yield correct results. 

Making use of the Wigner functions obtained in Paper I, in section ITT] we compute a CP-violating 
collision source in the scalar kinetic equation. The source originates from the fermionic Wigner 
functions, which we model by the CP-violating equilibrium solutions accurate to the first order in 
h. The source appears in the fluid continuity equation. We also briefly discuss the case of mixing 
fermions. This represents a first calculation of a CP-violating source in the scalar kinetic equation, 
and can be used, for example, for baryogenesis scenarios mediated by scalar particles, e.g. Higgs 
particles or squarks. 

CP-violating collisional sources in the fermionic transport equations are the subject of section 12.21 
We calculate the leading order CP-violating source, which appears first when the self-energies are 
approximated at the one-loop order. The source is of first order in gradients, it is suppressed by 
the wall velocity, and its parametric dependence is qualitatively different from the semiclassical force 



source. For example, the coUisional source is suppressed by the Yukawa couphng squared, and it 
exhibits a kinematic mass threshold - it vanishes when M < 2\m\, where M denotes the scalar mass. 
Moreover, the source appears in the fluid Euler equation, to be contrasted with the semiclassical 
force, which sources the continuity equation. The source we compute arises from a nonlocal (one- 
loop) contribution to the collision term, and it differs qualitatively from the spontaneous baryogenesis 
sources obtained so far in the literature y, |5|, ly, lD, |8|, |9J , where the source is computed from the (local) 
tree level mass insertions. 

In order to study equilibration, in section IHl we use a linear response approximation to compute 
relaxation rates for spin quasiparticles from the one-loop self-energies. When viewed in the Boltzmann 
collision term, these rates capture an on-shell emission or absorption of a scalar particle by a fermion, 
such that, as a consequence of the vertex energy- momentum conservation, the rates vanish in the 
limit of vanishing masses. Moreover, for finite fermion and scalar masses, the rates are kinematically 
suppressed, illustrating the shortcomings of the one loop approximation to the self-energies. Since 
the two-loop rates, which include the radiative vertex corrections to the above processes, but also the 
(off-shell) scalar exchange, do not vanish in the massless limit, they typically dominate relaxation to 
equilibrium, even when Yukawa couplings are quite small. To calculate these rates in the collision 
term of the Kadanoff-Baym equations would require an extensive calculational effort, and it is beyond 
the scope of this work. We begin section El by showing that thermalization rates do not depend on 
whether one uses equilibrium Wigner functions (approximated at leading order in gradients), or the 
standard spinors (with a definite spin). Next, we point out that in the ultrarelativistic limit the 
spin rates reduce to the (more standard) helicity rates, provided helicity is identified as the spin 
in the direction of propagation, which accords with our findings in section 11.51 This establishes a 
connection between the relaxation rates for spin and helicity states, which are almost exclusively 
quoted in literature. 

By taking moments of the transport equations obtained in section II. 5[ in section 0] we derive and 
analyze the fermionic fluid equations for CP-violating densities and velocities. We observe that 
the continuity equation contains the CP-violating source arising from the semiclassical force, while 
the Euler equation - which expresses the fluid momentum conservation - contains the coUisional 
source. The earlier works on baryogenesis, which were based on a heuristic WKB approach, have 
used helicity states, for which the semiclassical source appears in the Euler equation. Of course, 
helicity states can be obtained by a coherent superposition of spin states, such that, in principle, 
it should not matter which states one uses. In a quasiparticle picture, however, the information on 
quantum coherence is lost, and the dynamics of helicity states differ in general from that of spin 
states. Namely, at the order h, the interactions with a CP-violating Higgs background coherently 
mix the states of opposite helicity, which cannot be described in a quasiparticle picture. Formally, 
one can combine the Boltzmann equations for different spin and sign of momentum (in the direction 
of the wall propagation), such that the source in the fluid equations appears in the Euler equation, 
just as if helicity were conserved. These new equations correspond to taking different moments of 
the Boltzmann equations. Based on the above comments, it is clear that it would be incorrect to 
identify these as the fluid equations for quasiparticles of conserved helicity. In fact, this identification 
becomes legitimate only in the extreme relativistic limit. Taking different moments of the Boltzmann 
equations, such that the semiclassical source appears in different fluid equations, results formally in 
a priori equally accurate fluid equations. This observation lead us to conclude that these two sets of 



equations can be used as a testing ground for the accuracy of fluid equations, which could otherwise 
be tested only by a comparison with the exact solution of the Boltzmann equation. We perform this 
test in section HI and find that the two sets of diffusion equations give results that differ at about 10 
percent level. Finally, we show how, starting with fluid equations, one arrives at a diffusion equation 
which contains both sources, and thus allows for a comparison of the strength of the semiclassical 
and collisional sources at the one-loop level. Our conclusion is that unless diffusion is very effective 
the semiclassical force dominates. A more detailed look at the sources and rates from the one-loop 
collision term analysis indicates a mild suppression of the collisional source with respect to the flow 
term source. We expect that this suppression persists at higher loops as well. However, to get a more 
quantitative estimate for the strength - and thus relevance - of the collisional sources, one would 
have to include the contributions from the two-loop self-energies, which is beyond the scope of this 
work. 



2. COLLISION TERM 

In Paper I we have considered tree-level interactions of particles in a plasma with spatially dependent 
scalar and pseudoscalar mass terms. In this section we extend our analysis to include CP-violating 
contributions from the collision terms of both scalars and fermions. Our study is mainly motivated by 
an ongoing dispute in literature concerning the question whether the dominant CP-violating sources 
for electroweak baryogenesis arise from the flow term or from the collision term. While here we 
present just a first attempt to address this tantalizing question in a controlled manner, we beheve 
that our analysis, which is based on the collision term with the self-energies truncated at one loop, 
serves as a good indicator for the typical strength of the coUisional sources, and may serve as a good 
guideline for a comparison of the baryogenesis sources from the fiow term (semiclassical force) and 
coUisional (spontaneous) sources. Sorne, but not all, of the results discussed in this section have 
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already been presented at conferences [10|, lU , 

The simple picture of coUisional baryogenesis is based on the idea of spontaneous baryogenesis (l3l |. 
According to this idea, since particles and the corresponding CP-conjugate antiparticles with opposite 
spin see different energy-momentum dispersion relations in the presence of a CP-violating bubble 
wall, rescatterings have a tendency to preferably populate one type of particle states compared to 
the other, and thus source a CP-violating current. This current then biases sphalerons, leading to 
baryon production. This idea has been further developed and adapted to electroweak barvpeenesis 
mediated by supersvmmetric particles, mostly charginos and neutraUnos of the MSSM |^,y,la, i^ JSl 
and NMSSM |l4| . This approach is to be contrasted with the semiclassical force mechanism 1^ |l£ , 
E3, H H HTill discussed in detail in Paper I. 

In the classical Boltzmann equation the fiow term, which is first order in gradients, is balanced by 
the collision term of zeroth order in gradients, and hence both can formally be considered to be of 
the order hP in gradient expansion. We shall assume that our truncation of the fiow term at second 
order in gradients (formally at first order in h) can be considered equivalent to a truncation of the 
collision term at first order in gradients, so we will consistently neglect all higher order terms here. 

All Wigner functions in our work are on-shell and can eventually be reduced to particle distribution 
functions ()1.61 11.12|) . These distribution functions can be split up into a thermal distribution, the 
Bose-Einstein distribution nf in the scalar case and the Fermi-Dirac distribution rieq for the fermions, 
and a correction function 6n^ and 6n, respectively. Let us for the study of the collision term assume 
that these deviations from the thermal distributions are suppressed by at least one order in gradients 
in comparison to them. After expanding the sums nf + 6n'^ and n^q + 6n in the collision term, we can 
organize the contributions into two groups: those which contain correction functions, and those which 
only contain the equilibrium distributions. The terms in the first group are the relaxation terms one 
usually finds on the right hand side of a Boltzmann equation. They try to drive the deviations to 
zero and vanish for Sn'^ = 6n = 0. We will study them in more detail in sectional The terms of the 
second group owe their presence to the fermionic Wigner function, which contains first order effects 
besides the non-thermal distributions. These terms are present even for 6n'^ = 6n = and therefore 
precisely constitute the above mentioned coUisional sources for the correction functions. There is a 
caveat, however: by assuming that the deviations are suppressed by at least one gradient compared 
to the thermal distributions we ignore possible contributions from the dynamical CP-even deviations 



FIG. 1: The two-loop diagram corresponding to the Yukawa interaction H2.1() contributing to the 2PI effective 
action (|^ . 

from equilibrium, which are formally of the order 0{vwh^), and hence may also contribute to the 
source in the collision term (recall that in section 11.5.31 we discussed how a CP-even deviation from 
equilibrium can source a CP-violating deviation in the flow term ()I.5.136|) ). 

We shall work in the model of chiral fermions ipL,R ()1.1|) - which couple via Yukawa interactions to 
scalars 0, 

Ant = -^Ly(l)tpR-i'Riy<pyi'L 

= -^P {Pr ®y(t) + PL® (#)0 ^ , (2.1) 

where Pl,r = (1 -F 7^)/2 denote the chiral projectors, y is the Yukawa coupling, which is a matrix 
(vector) in the scalar (fermionic) flavor space, and ® denotes a direct product of the spinor and 
flavor spaces. Provided the momenta are not too small, the dynamics of fermions and scalars at 
the electroweak scale is expected to be perturbative, and the simplest nontrivial contribution to the 
effective action r2 in p.2.21|) from the two-particle irreducible vacuum graphs arises from the two- 
loop diagram shown in figure H A perturbative treatment is justified as long as we are interested 
in modelling the dynamics of nearly thermal quasiparticles, which are predominantly responsible for 
the sources and thermalisation processes discussed in this work, since a majority of quasiparticles 
carry nearly thermal momenta. After an easy calculation we find 

Ps = - /" (1% d% Tr [Pr ® y'S{u, v)Pl ® y''^ S{v, u)] Aw (n, v) , (2.2) 

from which we get the following one-loop expressions for the self-energies: 

5T2[A,S] 



I[1i';{u,v) = iac 



T,'"'(u,v) = —iac 



5A-(t;,«) 
^Tr [Pr ® y^'S^'^iv, u)Pl ® y'^ S''%u, v)] (2.3) 

sr2[A,s] 






= i[Pl® y'^S'^^u, v)Pr ® y''A'if\{v, u) + Pr ® y'S^\u, v)Pl ® /'Af,?(n, v)] . (2.4) 

The indices / and /' denote scalar flavor. For the treatment of the collision term we need the self- 
energies 

n<'>(M, v) = n+-'-+(M, v) = -iTi [Pr (g) yS>'<{v, u)Pl ® y^S<'>{u, v)] (2.5) 

S<'>(m, v) = E+-'-+(m, v) =i[PL0 y^S<'>{u, v)Pr ® yA>'<{v, u) (2.6) 

+ Pr0 yS<'>{u, v)Pl ® y^A<'>{u, v) 

8 



Equations (jl.!2H1.3|) together with (J2.3II2.4|) represent a closed set of equations for the quantum dy- 
namics of scalar and fermionic fields described by the Lagrangean p. HI . This description is accurate 
in the weak coupling limit, in which the Yukawa couplings satisfy \\y\\ ^ 1. One can in fact study 
the quantum dynamics of scalar fields also in the strongly coupled regime by making use of the 



1/N expansion, where N denotes the number of scalar field components [2j, |23|, |2^. Rather than 
studying this rather general problem, we proceed to develop approximations which are appropriate 
in the presence of slowly varying scalar field condensates close to thermal equilibrium. In this case a 
gradient expansion can be applied to simplify the problem. 

By making use of ()I.2.8j) it is easy to show that the self-energies 11* and 11* satisfy p.2.37j) with 11^^ = 0, 
representing a nontrivial consistency check of our loop approximation scheme for the self-energies. 
Of course this check works as well for the fermionic self-energies. 

2.1. Scalar collision term 

We begin with the analysis of the collision term ()1.4j) of the scalar kinetic equation in the Wigner 
representation. We expand the diamond operator to first order in gradients and obtain two contri- 
butions, in the following referred to as zeroth and first order collision term, respectively (although 
the "zeroth order" contribution in fact is of first order in gradients, as we will see): 

= -- (iu>iA< - m<2A>) + - o ( {zn>} {tA<} - {m<} {iA>} ) . (2.7) 

We shall now consider the scalar collision term for a single scalar and fermionic field, and then 
generalise our analysis to the case of scalar and fermionic mixing in the next section. 

2.1.1. The scalar collision term with no mixing 

The one-loop expression for the single-field scalar self-energy (vacuum polarization) ()2.5j) in the 
Wigner representation reads 

ffl<'>(fc,x) = -\y\' j "^'^^^'f {2nY6\k + k' - k")TT [PmS>'<{k' ,x)PLzS<'>{k" ,x)] , (2.8) 

where Pl,r = (1 =F 7^)/2 are the projectors on left- and right-handed chirality states. A graphical 
representation of this one-loop self-energy is shown in figure El The hermiticity property p.2.6j) of 
the fermionic Wigner function immediately implies the hermiticity relation 

{m<'>{k,x)y = iU<'>{k,x), (2.9) 

which holds in general for reasonable truncations of the scalar self-energy. With the above expression 
for the self-energy we can write the zeroth order collision term in ()2.7j) as 

C^o = % r-i^{2i.Y5\k + k'-k") 



2 J (27r)8 

Tr [PRiS<{k')PLiS>{k")\iA<{k) - Tr [PRiS>{k')PLiS<{k")\iA>{k)\, (2.10) 



■■^■o--^-- 



FIG. 2: The one-loop scalar vacuum polarization diagram H2.8|) from the Yukawa interaction H2.1() . which 
contributes to the collision term of the scalar kinetic equation H1.2|) . 



where we suppressed the (local) x-coordinate dependence for notational simplicity. We insert the 
spin-diagonal decomposition (jl.llj) of the fermionic Wigner function and take the trace. Because of 
the chiral projectors within the trace, the component functions gf and (^1 do not contribute. With 
the help of the relation 



Tr 



PRPAk')Ps"ik") 



^ ' 1 + s's" ° ° .^ J 

KqKq 



(2.11) 



we then find 



C, 



4 



E 



d^k'd^k" 



(27r) 



(27r) V(A; + A;' - A;")sV' l + s's 



/ // U U I 



y h" 



X { (^^'<(fc') - gi^k')) {g^;>{k") + gi'>{k")y/\<{k) 
- [gi'-ik') - gi>{k')) [gi'<{k") + gi'<{k")y/\>{k) 



(2.12) 



In our tree-level analysis of fermions in section 11.51 we found relations that enabled us to express all 
fermionic component functions in terms of g^. Here we need only equation p.5.67|) for g^, which for 
one fermionic particle reduces to 



^3 



s<,> 



1 

ko 



1 s<,> I ^ I |2/)/o s<,> I ^^ 

sk.gQ +-^\M Odk.go + -C^ss 
2ko ^ 



(2.13) 



In the following we neglect the collisional contribution in this relation, since it would result in terms 
which are of higher order in the coupling constant. We explained in the beginning of this section 
that the CP-violating sources, which contribute at first order in gradients to the collision term, are 
obtained by replacing the distribution functions in the on-shell Unsafe of the Wigner functions by the 
respective thermal distributions. For the remaining component ^^q of the fermionic Wigner function 
we therefore use 



5'oeq 



27r5 k^ 



1 



\7n\ 



ko 



^^>q = -27r5(k^-\m\^ + ^s\m\'e')\ko\{l-n,^Cko)) 

^ fcn ^ 



with the Fermi-Dirac distribution 



neq{ko 



1 



(2.14) 



(2.15) 
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instead of the full distribution function. Note that we are working in the wall frame, which moves 
with the velocity v^ within the rest frame of the plasma. Boosting the fermionic equilibrium Wigner 
function from the plasma frame to the wall frame results in the appearance of 



h = Jw{ko + Vyjk^) 



(2.16) 



in the argument of the equilibrium distribution. The CP- violating first order effects in the fermionic 
Wigner functions that have the potential to create coUisional sources are the derivative term in the 
relation ()2.13j) between g^ and Qq and the energy shift in ()2.14|1 . Since in the scalar Wigner function 
there is no first order correction other than the non-thermal distribution, we can here simply use the 
equilibrium form 



zA< (A;) = 2n5ik' - M2)sign(A;o)<(fco) 
tA>(k) = 27r6{e-M')sign{ko){l + ntJh)) 



(2.17) 
(2.18) 



again with the Bose-Einstein distribution 



ntik 



:./3fc0 



(2.19) 



evaluated at ko. All terms in ()2.12|) without derivatives acting on g^ then drop out because of the 
KMS relations (c/. section IT.2.5|) 



9oeq 



^l3ko A < 



eq 



and the energy conserving delta function 
12 r d^k'd^k" 



C,o = ^|m|V Y. 



X 



(27r)8 
s"k" 



{2iTy6\k + k' - k") s's" + 
s'k' 



y h" — h' ■ h" 
h' h" 



(2.20) 



(2.21) 



-^d,,y^<ik')(i - '-^)gQik") - (i + '-^)g'o^,ik')^d,,gQik" 



-^d,,9Z(k')(l-'-^)9Ctik''' 



L/tQ 



Jc" 



y 

l+'-^)9-^^,ik')±d,,gf<ik 



^Kik) 



fcn 



h" 



^Kik) 



This expression for the scalar collision term is already explicitly of first order in gradients, so we 
can neglect the spin-dependent shift in the spectral delta functions in ()2.14j) . The resulting spin- 
independent function, living on the unperturbed mass shell 6{k^ — |mp), is denoted by ^^oeq- We can 
now sum over spins and, taking into account correction terms caused by the momentum derivatives 

(2.22) 

use the KMS relation for the remaining terms to obtain 

d^k'd^k" 



dk.9oUk) = -e^"' {d,^g',<^ + P^^v^g'o^^) 
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(f>0 
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(27r)8 
9^oq{k')9lq{kyA<^{k) 



{2Ttf5\k + k' - k") 



1 



1 



(fco)^ 



(^oO 



(2.23) 
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This coUisional source is proportional to the wall velocity Vyj. This displays the fact that for a wall 
at rest the leading order Wigner functions indeed form an equilibrium solution. 
At linear order in v^, the equilibrium relation 



fi'Oeq ( ^o) ~ fi'Oeq (^o) 



(2.24) 



allows one to make the replacement k' *-> —k" in the last term within brackets, so that Eq. ()2.2Hj) 
can be rewritten as 



C<i>Q{k) 



^i^{k) 



\y 






T 

d^k'd^k" 



5\k + k' - k") 



(K 



1 h' h" — y ■ h" 

y h" 



9oetik')9&k"). 



(2.25) 
(2.26) 



tQn.Q 



Recall also that as a consequence of the KMS condition X?"q(/c)zA<(A;) = T^{k)iA^{k). In Ap- 
pendix El we show how, by performing all but one of the integrations in ()2.26jl . X<^o = ^m '^^^ ^^ 
reduced to a single integral ()A15|) . such that the scalar collisional source simplifies to 



a 



4>o 



-\y 

M2 



r,vJmr9'6(M-2\m\) f°° dk',,^, ,2 ^9, ,,, ,,,, , ^, , ,,,, 

'^^ ^ ^ / —9{ek'^ - /?2) e{ko)6{ko - u^) + e{~ko)Siko + u^) 

61 LO^ Jq UJq 



^- 



1 



+ 53- 



fc. 



-2|"^l ^/{k,-c^){c+-k,) ^{k,-c.){c+-k,) 



k' 



The c± and /3o are defined in Eqs. ()A14|) and ()A9|) in Appendix |X[ and 

^ , . 1 „ ^ , . 1 



Jq 



%q[^4 



^P^g> _ I 



/o = /o = n^qiuJo) 



QfSi^o + 1 ■ 



Note that the source contains a mass threshold: it is nonzero only if 

M > 2\m\ . 



(2.27) 



(2.28) 



(2.29) 



The physical reason for this threshold is that in the one loop self-energy only on-shell absorption and 
emission processes are allowed. Together with the energy-momentum conservation, this leads to the 
above condition. 

We now make a slight digression and consider the first order collision term in ()2.7|) . Because of the 
derivative in the diamond operator this term is already of first order in gradients and we can use 
the leading order equilibrium Wigner functions for its computation. Then the KMS relations ()2.2()j) 
induce a similar relation for the scalar self-energy. 



Il>(k)=e^'m<{k) 



eq\ 



(2.30) 



and we find 



C, 



</>! 



■-(d,Jko)d,(ll>^{k)zA<^{k) 



-^PlwVwd: 



vr o a I d'^k'd'^k ,^_-^i 



(2vr; 



3 {2T,)%k + k'- k")TT[PnS<^{k')PLS:^{k")]zA<^ik) . (2.31) 
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Following a similar procedure as in the leading order term, this can be reduced to 

-{2ny6{k + k'- k" f '''',': ■' 9Uk')9UkyK(f^) ■ (2.32) 



„,|2 r rl'^h'rl'^h" h' h" — h' ■ h" 



Ccj,i = -i—(3^wVwd, 



{2n) 



h' h" 



It is now clear from Eqs. (jl.9|) and p.4.36|) that this first order term contributes only to the constraint 
equation, since {C^i}'^ = —C^i is antihermitean. Furthermore, one can show that this contribution to 
the constraint equation results in a subdominant contribution in the kinetic equation, and hence can 
be consistently neglected. 

In anticipation of the discussion of fluid equations in section |3] we remark that the CP- violating 
sources in the fluid equations are obtained by taking moments of the scalar kinetic equation (jl.9|) for 
the CP-violating distribution function. We note that the zeroth order collisional source C(^o f!2.23p is 
real and symmetric in the momentum argument: 






(2.33) 



Qo = ^<t>o , C^o{-k, x) = C^o{k, x) . 

The symmetry is a simple consequence of the relation 

Kko + K- k'^)gl^{k',)g<^^{k'^) -. 6{-ko + k', - k'^) g<^^{k'^) g>^^{k',) , (2.34) 

where we made use of Eq. ()2.24j) . We furthermore recall that the first order source ()2.32|) is antiher- 
mitean, such that the following relation holds: 

^ [C^{k, x) + 4(A;, x) + C;{-k, x) + Cli-k, x)] = 2C^,{k, x) . (2.35) 

From Eq. ()A5j) and ()A8j) it follows that the integrand function in the scalar collisional source C^ = C^o 
is of the form j dk'^F{kzk'^, ^li^'z )•> ^^ ^^sX simple changes of variables imply 

dk, / dk'^F{kX,klk'^') = dkj dk',F{kX,kl,k'^') 

oo J Jo J 

/O p /"CxD p 

dk,k, dk'^F{kX,klk'J^) = - dk,k, dk'^F{kX,klk'J^), (2.36) 

oo J Jo J 

and we immediately conclude that the source in the scalar Euler fluid equation, which is obtained 
by taking the first moment of Eq. p.9|) . vanishes: 

(2.37) 






^4, 

The same conclusion is reached upon observing that C^ is symmetric under k^ — * —k^. The source in 
the continuity equation, which is the 0th moment of p.9|l . does not vanish, however. Indeed, upon 
performing the fco, k^ and angular (azimuthal) integrations in ()2.27|1 . we obtain 
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d'^k 

(2^ 



C. 



00 



\y\ 



iw' 



Vu,{\m\'e')Te{M-2\m\)l,^ 



-Lr-i 



Jo ^<t> Jo ^0 
M2 



X 



/ k' 

— — - arctan - — : 
2\m\ \\m\ 



Po 
k' 



\m\ ( k 

1 ; — arctan - — : 

k' \\m\ 



k' 



X /o*(l-/^) ^(c.^')(l-/o)k=.i-.,-^K)(l-/o')k= 



(2.38) 
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FIG. 3: The integral in the CP-violating colhsional source H2.38() of the scalar continuity fluid equation as a 
function of the mass parameters, M/T and \m\/T. 

We now evaluate this collisional source numerically as a function of the mass parameters, M/T 
and \m\/T (M > 2|?7i|). The result is shown in figure El The source peaks at \m\ k. IT and 
M ^ (3 — 4)T, and it is (exponentially) Boltzmann suppressed when |m|,M ^ T. Note that the 
nature of the collisional source ()2.38|1 is quite different from the source p.5.135|l arising from the 
semiclassical force in the flow term. In particular, the source fl2.38|l is a consequence of CP- violation 
in the fermionic Wigner function, which circles in the one-loop self-energy shown in figure El and it 
is hence suppressed by the Yukawa coupling squared, \y\'^. Moreover, the source exhibits a kinematic 
mass threshold, M > 2\m\, which is a consequence of the on-shell approximation in the fermionic 
Wigner functions. The threshold is of course absent in the semiclassical force source. To get a rough 
estimate of how strong the source can be, we take Vw ~ 0.03, \m\ r^ T, M r^ 3T, \y\ ~ 1, 6*' ~ O.IT 
(maximum CP- violation), such that S^"^^ ~ lO^^T^. 

2.1.2. The scalar collision term with mixing 

We shall now study the scalar collision term in the case of scalar and fermionic mixing, for which the 
Lagrangean ()2.ip reads 



£Vu = -0'^L2/Vii + h-c. = -V (Pr ® <\>'v' + Pl0 <p'*y'^) ^ , 



(2.39) 



where y = {y\j) is a matrix in the fermionic flavor space and / denotes the different scalar particles. 
The one-loop scalar self-energy is a straightforward generalization of ()2.8p . 

2n<;>(fc,x) 



{2'nY5^{k + k' ~ k")TT (Pn ® yHS>'< {k', x)Pl ® y^'US<^> {k", x)) (2.40) 

where the trace is now to be taken both in spinor and in fermionic flavor space. The mixing in the 
bosonic sector was handled in the flow term by a unitary rotation into the basis in which the mass 
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is diagonal (c/. section IT^. resulting in Eq. (|1.9p . so that an identical rotation has to be performed 
also in the collision term ()2.7j) . 

Let us consider first the influence of the scalar mixing. The scalar collision term is now to be 
understood as a product of two matrices in the scalar flavor space. After the rotation into the 
diagonal basis the form is maintained. The self-energy gets rotated as iH^'^ — > ^11^''* = UiH^'^W , 
which corresponds to the following redefinition of the couplings, y ^ yd = Uy and y"^ -^ yl = y'^U^ 
in ()2.40|) . The diagonal contributions to the zeroth order collision term contain the terms y\y\ iA^'f'i, 
such that the off-diagonals (1 — (5i;/)iA^j^ may give rise to a CP-violating contribution when either 
the Yukawa couplings |/^ or iA^^^J are CP-violating. Treating the latter contributions properly would 
require the dynamical treatment of the off-diagonals, which is beyond the scope of this work. Note 
that there is no contribution at order h coming from the diagonals of iA^'^ . For a diagrammatic 
illustration of the problem see figure ID below. The diagonal elements of the first order collision term 
are imaginary, and therefore do not contribute to the kinetic equation. 

Next, we consider the question of fermionic mixing (c/. section II.5.2.1|) . In this case the mass m 
in the free part of the fermionic Lagrangean (ll.lj) is a non-hermitean matrix in flavor space. The 
diagonal basis is defined by 

TUd = UmV'' = diag |mj|e*^-' , (2.41) 

where U and V are the unitary matrices that diagonalize m. Since the self-energy fl2.4U|) is a trace in 
the fermionic flavor space, its form stays invariant and we just have to replace the fermionic Wigner 
function and the coupling matrix by the rotated versions. For notational simplicity we shall work 
with only one scalar particle, and drop the superfluous indices d. The fermionic Wigner function 
may contain both diagonal and off-diagonal elements in the flavor diagonal basis. 
Let us first consider the contribution from the diagonal elements. To this purpose we write the 
fermionic Wigner functions as S^j'^ -^ 5ijSf'^ + 0(/i), so that we have: 

The trace is now to be taken in the spinor space only. We can go through the same steps as in 
the non-mixing case of section I2.1.H always keeping track of the indices attached to the fermionic 
functions. One difference in comparison to the non-mixing case is that in relation p.5.67j) between 
gl and f^g there is an additional term from the fermionic flavor rotation. The expression analogous 
to (IT^ is then 



C^o{k,x) = --f3v^J -^^-^{2nr6\k + k' - k") ".^.J' " (2.42) 



X 



V- , . , . „ . / (\m\\e' + 2A,)). (|m|2(^' + 2A,)), 



yi2 



where A^ = ^{VdzV^ — UdzU"^). Upon performing a similar change of variables as in the non-mixing 
case, k'^ <-^ —k'^ and i ^^ j, we have 

6\ko + k',- kl)gi>^{k',)gi'<^{kl) -. 5\-k, + k',~ A;(,')4<q(fc;')C> (^o), (2-43) 
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such that Eq. (J2.42J) can be recast as 



-00 



^M±M!l^«ij>_^,(,).,^<(,) 
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^3 



^^Oiii^) 



d^k'd^k" KK - K ■ K 1 

,, (27r)^^(fc + k'- k") ° V Jl L^>.'<(A:')(?o1'c,('^")-^ 



(2.44) 



This integral can be evaluated in part; upon performing analogous steps as in the one field case, one 
arrives at the result 



1^0ij{k) 



9{M - \mi\ - \mj\) 



k'dk'e{k^k'^ - (3l 



-. dk' k'' + kx - ^-^^-i-f -i-^-i^ 






fci^ \^> 



OiJ 



{K-K){K~K) 



X 



^K, - uj^) [e{-h)Ui - /^) + e{ko){i - /^j^y 
-^(o;,^ - J,,) [e{-k,)fJo, + ^(^o)(i - /;)(! - /o;-)] 



(2.45) 



which, in the one field limit, reduces to Eq. ()A10|) . By making use of the integrals ()A11|) . the k'^- 
integral can be evaluated. The final k' and k integrals have to be performed numerically, the results 
are very similar to the ones presented in figure El 

Like in the non-mixing case the first order collision term is imaginary and therefore does not contribute 
to the collisional sources at order h. The self-energy stays imaginary, and the derivative from the 
diamond operator acts on the self-energy as a whole and does not notice the rotation. Finally, we 
remark on the contributions from the off-diagonal elements of the Wigner functions, (1 — 5ij)iSf-'^ ^ 
which are suppressed by at least one power of h. In the case when the fermionic interaction and fiavor 
bases do not coincide, it is clear that the off-diagonal elements of iS^'^ could in principle contribute 
to a CP- violating source at order h in the bosonic collision term ()2.7j) . A proper treatment of this 
problem would require a dynamical treatment of the off-diagonal elements of ^5'^^''* (in the flavor 
diagonal basis) and it is beyond the scope of this work. 



2.1.3. Discussion of scalar collisional sources 

Contributions to the scalar CP-violating collisional source could, in principle, arise from the scalar 
or from the fermionic sector of the theory (jl. 112.11) . In this section we have shown that, to order h in 
an expansion in gradients, there are no CP-violating contributions to the collision term of the scalar 
kinetic equation arising from the scalar sector. Our proof includes the more general case of mixing 
scalars. A simple qualitative argument in support of this is shown in flgure lU 

On the other hand, as we have argued in sections 12.1.11 and 12.1.21 the fermionic contributions to 
the scalar one-loop self-energy do result in CP-violating effects which contribute to the collision 
term of the scalar kinetic equation. As illustrated in flgure El these effects arise from CP-violating 
contributions to the fermionic propagators which run in the self-energy loop, and can be thought of as 
radiative contributions to the CP-violating scatterings off the Higgs condensate. These contributions 
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0=''0- 0--0 



FIG. 4: The one-loop collision term in the scalar kinetic equation in a theory with the Yukawa interac- 
tion H2.1() . The interactions of the scalar propagator with the scalar condensate are of the same nature as 
the tree- level interactions discussed in section 11.41 Just like in the tree-level case, at order h, scalars do not 
feel the condensate. 





■■■<->* 



FIG. 5: The collision term in the scalar kinetic equation as in figure |3 The interactions of the fermionic 
(loop) propagators with the scalar field condensate that can contribute in CP-odd manner at the order h to 
the scalar kinetic equation. 



can source CP-violating scalar currents in the Higgs or sfermionic sectors, for example, and thus may 
be of relevance for electroweak baryogenesis. From our one-loop study of the CP-violating source 
shown in figure 121 it follows that: (a) the source has a very different parametric dependence than the 
semiclassical source shown in figure 11.51 and (b) the source tends to be large in the case of fermions 
whose mass is of the order of the temperature, and where the scalar mass is a few times bigger. With 
T ~ 100 GeV, popular extensions of the Standard Model, which include supersymmetric theories, 
typically contain such particles. 

2.2. Fermionic collision term 

Now we turn our attention to the collision term for fermions. Again, we expand in gradients, and 
distinguish between the zeroth and first order collision terms. 



where 



C 



ipO 



C-ip — C^o + C^i + • • 



S>5<-S<5> 



C^i = -^o({S>}{S<}-{S<}{S>}). 



(2.46) 

(2.47) 
(2.48) 



CP-violating contributions to the collision term appear, like in the scalar case, at first order in 
gradients. Since the calculational procedure of the sources is very similar to that in the scalar case, 
we shall outline only its main steps. 
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s ; \ s 

FIG. 6: The one-loop fermionic vacuum polarization diagram (c/. Eq. H2.49() ') from the Yukawa interac- 
tion ()2.1|) , which contributes to the collision term H2.46() of the fermionic kinetic equation (|1.3|) . 

2.2.1. Fermionic collision term with no mixing 

The single- field fermionic self-energy, when approximated by the one-loop expression ()2.4|) and written 
in Wigner space, reads 

X [i^>'<{k")PLtS<^>{k')PR + i^<'>{-k")PRiS'^'>{k')PL]. (2.49) 

A graphical representation of this self-energy is shown in figure IHl We again use the spin diagonal 
Ansatz fll.lll) for the fermionic Wigner function, which in the self-energy shows up between two chiral 
projectors. This can be simplified by 

PLY.Sr{k')Pn = Y.'PAk')s'P,[^'r'%^^^{k')-^^gi<>{k')\ , (2.50) 

s' s' 

and similarly for Pl ^^ Pr, such that the self-energy becomes 

(2.51) 
Again we use thermal distributions for the Wigner functions, and with 

zA< (-A:") = ^A> (fc") (2.52) 

the chiral projectors drop out of the zeroth order collision term: 

C^o = ^Ey -^^5\k-k' + k")J2PAk')Ps{k)s' (2.53) 

X [tA<^{k") (7V^1> (^') - I'gZik')) {s^'l'g^Uf") - ^^'&W + ^ieq(^) - ^l'92Uf') 



It is not hard to see that, at leading order in gradients, the collision term vanishes. 
The further procedure is now analogous to the scalar case. We use equations p.5.65III.5.67j) . simplified 
to the case of only one fermionic particle, to express the functions g^^^ (i = 1, 2, 3) in terms of ^^oeq- 
The terms without derivatives acting on g^ vanish because of the KMS relations ()2.20|) and the 
energy-momentum conserving (5-function, while Eq. ()2.22p shows that all surviving terms are wall 
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velocity suppressed - as they should be: 
12 r d^k'd^k" 
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(2.54) 



Since this expression is explicitly of first order, we switched to the spin-independent functions ^foeq- 
Let us now have a look at the first order collision term ()2.48|) . The fermionic equilibrium self-energy 
satisfies the KMS relation 



^'^..{k) = -e^^°E< (fc) , fco = l^h + v^kz) 



so that - just as in the scalar case - we can write: 



C^i = --p'jtoVutdz (S<S'>) . 



(2.55) 



(2.56) 



Eventually we are interested in the contribution from the collision term to the fermionic kinetic 
equation p.l9|) . To obtain this contribution we have to multiply the collision term with the spin 
projector Ps{k) and take the trace. We find 



Ti [PAk')Psik)] =l + ss^ 



koK 



(2.57) 



while the traces containing additional factors of 7^ and 7^ vanish. 

The trace of the first order collision term is imaginary and so does not contribute to the kinetic 

equation. We find 



mk) = -mT[iPs{k)c^] 

\y\ n I \2n' 



d^k'd^k" 
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(2.58) 



This collisional source is symmetric under k ^ —k, therefore it contributes to equation ()1.19j) for 
Sf^, where the collisional contribution is 
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(2.59) 
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but not to equation (J1.18J) for 6f^. To further simplify this expression we keep only the terms that 
contribute at linear order in the wall velocity Vw The (5-functions in goeq and iA^q project onto the 
classical shell, ko = icuo, K = ^^o ^^^ k'^ = ±uj'^. Together with the energy-momentum conservation 
(expressed through the delta function in Eq. ()2.58p ). this leads to the constraint 

loquj'q - Iml"^ > k ■ k' . (2.60) 

After multiplying by —2 and adding fc^ + k''^ this can be rewritten as 

(coo-cu'^f < {k-k'Y < {k-k'Y + M\ (2.61) 

where we assumed M^ > 0. With the spatial part of the momentum conserving ^-function we 
recognize the last term as the square of u'^ and finally conclude 

±{uo-iUo)<uj';. (2.62) 

This means, however, that after performing the feg; ^o ^^"^ ^" integrals the contributions with 6{u!o — 
u'q ± u'!) vanish, as well as the ones which contain 6{±{u!o + u!q + uj'I)), such that we find 



ozvr J ujquJqUJ^ 



X 



X 



Kko + ^o)/o/o(l + /oH - ^(^0 - uj,){l - /o)(l - f,)f, 
(^ + -k-^(^^^o + h-kd)- (2-63) 

\UJqUJq UJqUJq J 

With uj"^ = ijJq \ ijj'q, the following relation holds 

(l-/o)(l-/^)/o" = /o/^(l + /o1, (2.64) 

and after inserting this into ()2.63|) . we immediately see that the collisional source is 

^0 = -7;7^s|mp6''/3t;^[(5(A;o + ^o) -"^(^o-^o)] 

OZTT 

d'k' r^p-6{u;o + u;',- c.^)/o/^(l + /o^) k'^ + k. \,, " • (2.65) 

Just like in the scalar case, parts of the ^'-integration can be done analytically. The calculation is 
performed in Appendix O 

In section 0] we derive fluid equations for the CP-violating part of the fermionic distribution function. 
For this we need the zeroth and first fc^-moment of the collisional source. In appendix O we show 
that the collisional source is odd under kz ^->- k'^, so the zeroth moment vanishes. The first moment 

^'^ ^'-sS-\k) = 4 / 7^^^^ = v^s\m\'eV-J^ . (2.66) 



{^nfuJo ' ' y,„>o (2vr)4 ^0 "^ - ' ' - 167r3|m|2- 

is nonzero, however. In figure [7| we plot the dimensionless integral X as a function of the fermion 
mass |m|; the scalar mass is chosen to be proportional to \m\. The source vanishes for small values of 
the masses, suggesting that the expansion in gradients used here yields the dominant sources. Note 
that the source contributes only in the kinematically allowed region, M > 2\m\. When the masses 
are large, |m|,M ^ T, the source is, as expected, Boltzmann-suppressed. 
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2.2.2. Fermionic collision term with mixing 

For several mixing scalar and fermionic particles the procedure is very similar to what we have done 
in section l^.l.2l for the scalar collision term: we apply flavor rotations to diagonalize the mass matrix 
in the flow term of our equations and study the effects of these rotations on the collision term. The 
mixing Lagrangean ()2.39|) yields a fermionic self-energy of the form 



S<'>(/t,x) 



d^k'd^k" 



(zAfi;<{k",x)PL®y'^S<'>{k',x)Pn®y' 
+ tA^i>{k", x)Pr ® y''S<^>{k', x)Pl ® y'^ 



(2.67) 



The quantities with scalar flavor indices form a trace and so their form does not change under the 
scalar flavor rotation, we only have to replace y -^ yd and iA -^ iAd. We again neglect the off- 
diagonal components of the scalar Wigner function, iAd^i>i -^ SwiAd^i, and so the self-energy simply 
becomes a sum of self-energies without mixing, with one contribution for each scalar mass-eigenstate. 
Since scalar mixing only appears within the self-energy, the derivatives appearing in the first order 
collision term do not change this result. 
It is therefore sufficient to study fermionic mixing with only one scalar particle. After the fermionic 
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FIG. 7: First moment of the collisional source (|2.65|) as a function of the rescaled mass |m|/r. The scalar 
mass is set to be a multiple of |m|. 
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flavor rotation, the diagonal components of the zeroth order collision term read 

^^^i2nr5\k - k' + k") 



(2.68) 



^E (fe^^^^(^"'^)^^^f (^''^)^^ 



+ y*^y,,zA>{k", x) Pr S>{k', x) Pl] S<{k, x) 



Here y and S denote already the rotated quantities, and we neglected the flavor off-diagonal elements 
of the fermionic Wigner function, Sij -^ SijSi. With the same steps as in the non- mixing case we 
arrive at the expression analogous to ()2.58|) : 



/^Oii 



■^(^IwV. 



^ , -4A^{'^^Ynk-k' + k'')iA>{k'') 



27r 



(2.69) 



7 , 2 dojeqi^ )9oieq{k) 

J 

J2 3?TrP,,(fc')P.(A;) [ - .^1 [\m\'{9' + 2A,)]^ + s'^^ [\m\'{9' + 2A, 






kok' 



Like in the non-mixing case the contributions from the leading order functions gf and «A^ vanish, 
if they lead to S{u!oi — u'qj ± u'^) or S{±{ujQi + uj'q- + cj^)). We find 



/C; 



Oil = 7/^^t" ["^(^o + ^Oj) - Kh - ^Oi)] 



(2.70) 
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TT 



(27r)3 2 
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\y^3? + \yji? 1 



LO 



s[ \m\ 



\2fQl 



k' 



k^ ^oi^'oj + k\\ ■ fc|| 



+ 2A,))^ J , + s{\m\\e' + 2A,)) 



i^Oil^Oi^i 



Oj 



^ CJoi^Oi^Oj ^i 



Oj 



From here on one can proceed as we did in the non-mixing case and perform some of the integrals. 
This is shown in appendix O 

In the first order collision term for the fermionic equation the diagonalization leads to extra terms 
with spatial derivatives acting on the rotation matrices. But none of them contributes to the relevant 
trace for the kinetic equation for Qq. 
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3. RELAXATION TO EQUILIBRIUM 

In the beginning of the last chapter we explained that the collision term consists of two parts: the 
sources, which are caused by the structure of the fermionic Wigner function, and relaxation terms 
due to the deviation of the distribution functions from thermal distributions. This chapter deals 
with these relaxation terms. In a first section, considering as example the scalar collision term, we 
show that the relaxation part of the collision term is equivalent to the right hand side of a usual 
Boltzmann equation in our one-loop approximation. In a second section we give explicit expressions 
for the relaxation terms relevant for the kinetic equations ()1.9j) and (jl.lHflTmUl . and comment on 
their consistency with the flow terms. In this chapter we work with only one scalar and one fermionic 
particle. The extension to several mixing particle species is quite straightforward. 

3.1. Boltzmann collision term 

Since we only keep terms up to first order in gradients in the collision term, we can for the sake of 
calculating the relaxation term use the leading order expressions p.2.66flT.2.67|) for the scalar Wigner 
functions, however with the thermal distributions nf^ = l/(e'^'^° — 1) replaced by the function n*^ = 
nf + Sn^, which is a sum of the thermal distribution and an unknown first order correction, 6n'^ = 
6n'^{k,x). For the fermionic Wigner function we use the spin-diagonal expressions (c/. Eqs (11.2.711 - 

OH) 

S^{k,x) = 2Tciy^Ps{k){p + rrih — i'y^'ma)6{k'^ — \7n\'^)sign{ko)ns{k,x) (3.1) 

s 

S>{k,x) = -2m^Psik){p + mh-i-f^ma)5{k'^-\m\'^)sign{ko){l-ns{k,x)), (3.2) 

s 

where 

Us = —I h 5ns , ko = -fu,{ko + v^k^) (3.3) 

gPfco _|_ 1 

denotes the Fermi-Dirac distribution plus an unknown, spin dependent, correction, 6ns = Sns{k,x). 
Note that the spin projection operator per construction commutes with ^. 

Since the scalar first order collision term contains an explicit derivative, any first order correction 
from the distribution functions would lead to second order terms and therefore can be neglected. 
With the above Wigner functions the scalar zeroth order collision term ()2.im) becomes 

C,o = -^ / ^^|^(27r) V(A; + k' - k") J^ Tr [PRPAk') f P,\k")f 



s' 



'2 \^\2\-„„nJ\ n^xih.ii'^ 



27T6{k" - |m|")sign(A;^) 27T6{k"' - \m\')sign{k'^) 2TT6{k' - M')sign{k 



0) 



ns'ik'){l-ns"{k"))n'^{k)-{l-nsik'))ns"{k"){l + n'^{k)) . (3.4) 

We introduce the notation 

Fs^s'ik', k") = Tr [PnPsik') f Ps"{k")fj (3.5) 

for the trace appearing in this expression and note that this object is antisymmetric in both of its 
momentum arguments. Upon having a closer look at the Boltzmann equation p.9|) for scalars that 
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FIG. 8: The scattering diagrams for an incoming scalar particle. Particles and momenta run from the left 
to the right. 



we derived in section II.4.H we see that the coUision term appearing on the right hand side effectively 
is 

dkoC^{k,x) , (3.6) 



1 

TT./O 



where C^ is defined in ()2.7|) . After performing the fcg-and feg -integrations and making some rearrange- 
ments, the relaxation contribution to ()3.4H3.6|) can be written as 



n 



dkoC^o{k,x) 






(3.7) 



X [{2iTy6\k + k'- k") [fs'+{k'){l - fs"Ak")fHk) - (1 - fs'Ak'))fs"Ak")i^ + fHk)) 
+ (27r)V(A; - k' - k") [(1 - fs'-{k')){l - fMk"))f{k) - f,^{k')f,.S"){^ + fm 

+ (2vr)V(fc -k'+ k") \{i - f,^{k'))f,,^{k")f^{k) - fs'-{k'){i - fs"-ik"m + fm 



where fs± and f^ are the on-shell densities for fermionic particles and antiparticles, and scalar 
particles, respectively. 

Let us now, for comparison, consider the collision term one would usually expect for a scalar particle. 
In figure |H1 we show all tree level processes for an incoming scalar particle given by the interaction 
Lagrangean 

Cmt = -^ {PRy<f) + PLy*<l)*) ^ . (3.8) 

The matrix element for the first process, the transition of the scalar particle into a fermion-antifermion 
pair, is 

MS''^ = -%yu,,{}^')PRV,{}^), (3.9) 

with the notations and conventions for the spinors as used in Ref. 25|. Calculating the absolute 
square of the matrix element is done as usual, but with the complication that we do not sum over 
spins. We find 



l-^i'i"(^',m' = bl'Tr \PnPA^) If' PAk")f] = \y\^Fs>Ak\k"). 



(3.10) 



This does not look immediately like an absolute square, since by sending k' to —k' the sign of the 
trace changes, but we have to keep in mind that the particles have to be on-shell and the above 
expression has to be multiplied by the energy-momentum conserving 5-function, 5^{k — k' — k"), so 
that eventually ()3.10|) is positive definite (one can also check this by actually performing the trace). 
In a Boltzmann collision term this matrix element appears with a minus sign and is multiplied by the 
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distribution functions (1 — /+)(! — /-) for the outgoing fermion and antifermion, respectively, and f^ 
for the incoming scalar particle. We can identify exactly this contribution in the second line within 
the curly brackets in Eq. ()3.7|) . Subtracted from it is the inverted process where a fermion and an 
antifermion annihilate to form a scalar (c/. Ref. [26| for similar analyses). In the same way we can 
find the contribution for the process shown in|Hl^b), absorption of a scalar by a fermion, and|Hfc), 
absorption by an antifermion, in the first and third line of ()3.7p . respectively. Subtracted are the 
inverted processes, where the scalar comes out. The process shown in|Sl^d) is kinematically forbidden 
and therefore does not appear in ()3.7j) . So the relaxation part of the scalar collision term is exactly 
what one would expect to find on the right hand side of a classical Boltzmann equation for a scalar 
particle. It can be shown that the same is true for the fermionic collision term and the Boltzmann 
equation for a particle with definite spin. 

3.2. Relaxation rates 

In section ir.4. II we derived the kinetic equation ()1.9|1 for the scalar CP-violating particle density Sfip. 
Obviously, expression ()3.4j) for the relaxation term is real, so that the right hand side of the kinetic 
equation becomes 

C^{k)+C^{-k) . (3.11) 



We insert fl3.4|l and send k' — > —k' and k" -^ —k" in the second term. The trace stays invariant 
under this change, so that we find 

dkp 

271 



C^{k)+C^{-k) 
'^^^{2T^Y5\k + k'- k") J2 Fs'Ak', k")2 j^ dk 



2 

6{k'^ - \m\^)sign{k'Q) 6{k"" - |mnsign(fc;') 6{k^ - M^)sign{ko) 
Usik') (1 - nsn{k")) n'^ik) - (1 - Us'ik')) ns"{k") (l + n^{k)) 

ns'{~k') {l-ns"{-k"))n'''{-k) + {1 - ns'{-k'))ns"{-k") (l + n'^(-A;)) | . (3.12) 



The terms in the last line can be rewritten as antiparticle densities based on the relations (jl.Sj) 
and (|1.15|) . Then we split up the distribution functions into a thermal distribution plus a first order 
correction term, 

nHk) = <(A;o) + Sn^ik) , n'^'^^ik) = nt^{ko) + 5n'*''%k) (3.13) 

ns{k) = n,^{ko) + Sn,{k) , nf{k) = n,^{ko) + Snf{k) , (3.14) 

and linearize with respect to the correction terms. This constitutes a rather standard linear response 
approach to equilibration, which is believed to lead to quantitatively correct results for systems close 
to thermal equilibrium, provided they are not prone to (hydrodynamic turbulent) instabilities J27|. 
Note that the Bose-Einstein and Fermi-Dirac distributions in (|3.13riXT^ are evaluated at ko = 
1w{ko + Vujkz), since the wall frame moves with respect to the plasma. Under extensive use of the 
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explicit form of the equilibrium distributions (|2.15p . (|2.19|) and energy-momentum conservation we 
can rewrite the first term in curly brackets as 



(l-n,^CK)) <(M (l-^eq(A;;)) ' 6nAk' 



^eq(A)[,)nfq(^o)neq(A;o) Sn,„{k") 



+ rie^CK) (l - neciCK)) (l + <(^o)) ' Sn^ik) . 



(3.15) 

Exactly the same manipulations are then applied to the CP-conjugate densities, and since the equi- 
librium distributions are spin independent, we can combine the two expressions to 



dko 
~2^ 



C^{k)+C^{-k) 

d^k'd^k" 



(3.16) 



\yr / dko 

'0 



{2TiY5\k + k' - k")J2^s's"{k', k") 



(27r)6 

X 5{k'^ - ImH sign(A;[,) 5{k"^ - \m\^) sign(A;^') 5{e - M^) sign(A:o)neq(A;[,)(l - neq(fc^'))<(^o) 
5n''{k') - 6n-^'{k') 5n'" {k") - 5n-f {k") 5n'^{k) - 5n%{k) 



X 



+ 



neq(fco) (l - '^eq(fco)) "^^M) (l " ^eq(fco)) '^eq(A;o) (l + '^eq(fco)) 



For the fermionic collision term the same procedure can be applied. For the combination of collision 
terms we eventually need in the kinetic equation ()1.19|) for the CP-violating distribution function 
5/", we find 



dkn 



TT 



IClik) + lCl{-k) 



\y\ 



d^k'd^k" 



E 



dkn 



X 6{k'^ - \m\^) sign(A;[,) 6{k"^ - \m\^) sign(A;;') 6{k^ - M^) sign(A;o) 
X [{27:Y5\k -k' + k")F,Ak', k")n,Xko) U - n,^ik',)) ntM) 



(3.17) 



X 



+ {27,Y5\k -k'- k")F^s-Ak', k")n,^{ko) [I - n^^k'^)) (1 + <(A;^') 
6n'{k) - 6n-;{k) 5n'\k') - Sn~/ {k') 5n^{k") - 6nHk") 



+ 



neq(A;o) (l - nec^iko)) n^qik'^) (l - neq(A;^)) n|'q(A;^') (l + n|'q(A;[,')) 



The on-shell ^-functions in these two equations finally project the differences of the deviation functions 

6n — dricp onto the CP-violating distributions Sf^ and 5fg. Of course, the relaxation term of the 

kinetic equation for (5/J is 

dk V 1 

° icm+^o'{-k) 



IT 



(3.18) 



so that here we find Sf^ as well as on the right hand side. This shows that our kinetic equations ()1.9|) . 
fll.l8|) and p.l9p indeed form a closed set of equations, or in other words, that in the relaxation terms 
of the kinetic equations for the functions Sf^, 6fg and 5/" only these functions themselves do appear. 
It might be surprising that the relaxation term for the scalar equation contains the fermionic (5/", but 
not Sf^. The reason for this is the occurrence of the chiral projectors in the interaction Lagrangean. 
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The relaxation terms that we have discussed so far can be shghtly simphfied further by performing 
the sums over those spins that do not appear in the deviation functions. But essentially the form 
for the rates ()3.16|) and ()3.17|) is as simple as it gets, since we have to integrate over the unknown 
functions 6f. In the next section we make a fluid Ansatz for the CP-violating distribution functions 
in order to facilitate numerical solutions of the kinetic equations. The fluid Ansatz expresses the 
distribution functions in terms of two unknown functions, the chemical potential and the plasma 
velocity, which, in the stationary limit and in the wall frame, depend only on the spatial ^-coordinate. 
The momentum dependence is completely specified by the fluid Ansatz, so that the integrals can be 
performed and result in interaction rates, which multiply the chemical potentials and plasma velocities 
of the respective particles in the relaxation term. In appendix ID] we show the results for the rates 
for the fermionic quantities in the equation for 6fg. 

The rates we obtain from the one-loop self-energies do not always capture the physically dominant 
processes responsible for thermalization, however. This is so because we calculate them by forcing the 
outgoing particles on-shell, which, as a consequence of the energy-momentum conservation, results 
in a kinematic suppression. Because of this suppression, the rates vanish as any of the (tree level) 
masses approaches zero and, in general, they are suppressed when the masses are small, that is when 
\m\, M ^T. The processes described by the one loop approximation in figure |H1 constitute absorption 
and emission processes only. The one-loop approximation completely misses 2-to-2 particle elastic 
scatterings, which are essential for particle transport. 

When the self-energies are approximated at two-loop level, apart from radiative vertex corrections 
to the one-loop rates, one also captures (tree level) 2-to-2 scatterings rates, which involve off-shell 
(boson or fermion) exchange. These rates, even though suppressed by two more powers of the 
coupling constant, do not in general vanish when (any of the) masses approach zero, and hence often 
dominate thermalization rates in a linear response approximation. For details of thermalization rate 
calculations starting with an (on-shell) Boltzmann collision term, we refer to Refs. [la, ll^ |2^ U^ 
l8fll |. A complete two-loop treatment of the collision term in the Schwinger-Keldysh formalism is a 
much more complex task. Nevertheless, we can say something constructive on results of such an 
undertaking. Indeed, our findings in the first part of this section indicate that the relaxation rates 
calculated from two-loop self-energies (with distribution functions projected on-shell) can as well 
be obtained by writing down the naive Boltzmann collision term. Then the scattering amplitude 
in the collision term is calculated from the 2-to-2 scattering processes, with the important caveat 
that one must not average over spins (or helicities), as it is usually done in literature, but needs 
the expressions for the eigenstates of the spin operator ()1.5.18I1.10|) . (For an example of a helicity- 
flip rate calculation, which does not involve averaging over helicities, we refer to Appendix B of 
Ref. [,16]-) Recall however that in the super-relativistic limit the spin operator S^ approaches the 

helicity operator p.5.22j) (multiplied by k = k/\k\). The rates involving particles with a definite spin 
are therefore in this limit identical to the rates for particles with definite helicities, defined as the 
projection of spin onto the direction of particle motion. 
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4. FLUID EQUATIONS 

In Paper I and in the previous sections we have shown that in the semiclassical hmit the equations 
of motion for the scalar and fermionic Wigner functions can be reduced to on-shell conditions and 
Boltzmann equations. In the first part of this section we closely examine the physical interpretation 
of the quantities we have been dealing with. Then we examine the highly relativistic limit of our 
Boltzmann equations for the fermionic particles and compare with the corresponding equations for 
particles with definite helicity. In the third part we make a fluid Ansatz for the CP-violating dis- 
tribution functions. The resulting fluid equations are numerically tractable, and we solve them in a 
very simple model. Finally we reduce the fluid equations to diffusion equations in order to make an 
estimate of the importance of coUisional sources in comparison to the flow term sources. 

4.1. Currents 



We begin this section by establishing some links between the quantities we dealt with in the previous 
sections and physical, measurable quantities. The best way to do this is to study various currents. 
The expectation value of the scalar current operator 

(4.1) 



j^{x) = (i(j)\x) d (j){x) 
can be written in terms of the Wigner function: 

j^{x) = -{du - d^)A<{u,v) 



l2^ 



kiA<{k,x) 



(4.2) 



This current contains information about the density of scalar particles (minus antiparticles) at some 
space-time point x. With the on-shell Ansatz (ll.(i|l for the scalar Wigner function, the components 
of the current can be written as 



jS(^) 



J<p[x) 



d^k 
(2^ 
d^k k 



(27r)3u;^ 



ft{k,x)-ft{k,x)-l 
ft{k,x)-f^{k,x) 



d^k 

(2^ 



5f{k,x)-l 



(4.3) 
(4.4) 



This confirms that f^{k,x) indeed measures the density of scalar particles and antiparticles with 
momentum k. The one in the first equation is the vacuum contribution and is therefore ignored. 
Similarly, the expectation value of the fermionic current operator 



4(x) = {^Pix)ri^ix)) 



which can be rewritten as 



Ji,ix) 



-Ti[^^'iS<{x,x)] 



d^k 

(2^ 



TiYiS<{k,x) 



(4.5) 



(4.6) 



contains information about the density of fermions. A careful consideration of the zero- component 
and the vector part of this current, and the three-component of the spin-density 



h 

5,(x) = -(^(x)7V^(x)) 



(4.7) 
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expressed in terms of the on-shell distribution functions (|1.14j) and (|1.15|) tells that fs±{k, x) measures 
the density of particles and antiparticles with momentum k and spin s (in the direction specified 
by ^[M^)- 

When considering electroweak baryogenesis, the final step in the mechanism that produces the baryon 
asymmetry is the sphaleron process: it turns the CP-violating flows, whose kinetics we study in this 
paper, into a net baryon density. To be more precise, the sphaleron acts only on the left handed 
SU(2) doublets, which are eigenstates of the chirality projector Pl, but leaves right handed particles 
unaffected. The densities of left- and right-handed particles 

Jl/Ri^) = {^l/ri'^l/r) (4.8) 

are linear combinations of the vector density and the axial density j^ = (%l)'~^^'^^ipy. 

3l,K = \{fT3'). (4.9) 

Finally, the vector density and the axial density are related to the fermionic CP-violating densi- 
ties (J1.16p and (jl.l7|) we defined earher, 

^:;^3^fsik,x) (4.10) 



fix) 



E [4^s\'-^f:(k,x) + s^mk,x)] , (4.11) 



which explains the terminology. We omitted the vacuum contribution to the vector density. It 
might seem strange that the equilibrium distribution contributes to the axial density, but this is 
just a consequence of defining f^"^ as the projection of neq to the mass shell ujs- Had we defined 
f^^ = neq(co'o), then the axial density would be expressed by 6f^ alone, but in turn there would be 
additional source terms in the Boltzmann equation for this function. 

Both CP-violating densities contribute to the left-handed density, so both could lead to baryon 
production. Recall however that in the equation (jl.lSj) for (5/J there is no source, neither from the 
flow term nor from the collision term. We therefore assume that this density, and consequently the 
vector density, vanishes and then we have j° = —j^ = —j^/2. The source in the equation for 6fg 
is proportional to spin, so that any contribution in 6fg caused by this source is also proportional to 
spin. Therefore the spin summation in the expression for j^ does not cancel the contributions against 
each other, but adds them up, wherefore the source in the equation for Sf^ eventually is an effective 
source for electroweak baryogenesis. The source in the scalar equation can contribute to baryogenesis 
only indirectly. It first creates a CP-violating scalar density, which then may be transferred into the 
fermionic CP-violating distribution function via the collision term of the scalar kinetic equation. 
Finally, a word on mixing particles. The sphaleron acts on the particles in the electroweak interaction 
basis, while the densities we are dealing with are those of particles in the mass basis. So in principle 
we would have to rotate these densities into the interaction basis. This is not necessary, however: 
the sphaleron is effective only in the symmetric phase, where the Higgs vacuum expectation value 
vanishes, and in those regions of the bubble wall where the expectation value is small. But for 
a vanishing Higgs expectation value the mass basis and the interaction basis coincide. So we can 
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compute the contribution to the axial density for each particle, using the densities in the mass 
eigenbasis, and sum them up. 

In the following we will not consider the equilibrium contribution to the axial density any further, 
another consequence of the sphaleron's refusal to work anywhere else than in the symmetric phase. 
Electroweak baryogenesis is most effective when there are CP-violating flows that are efficiently 
transported away from the wall into the symmetric phase, where they are converted into a baryon 
asymmetry. The equilibrium contribution to the axial density, however, sits right on top of the wall 
and does not participate in diffusion. 

4.2. Spin or helicity? 

In the highly relativistic limit we can approximate (Do by \kz\. In this case the density of left-handed 
particles can be written as 

where we allowed for a possible contribution from the vector density. But this is of course precisely 
the density of particles with negative helicity minus the density of antiparticles with positive helicity 
(for massless antiparticles chirality is opposite to helicity): 

j° oc / dh (5/-+ - 5/+_) + / dh (5/++ - 5/__) . (4.13) 

Jo J -oo 

In fact, in every earlier work on EWB that used the WKB method, the Boltzmann equations were 
written directly for the densities of particles with a definite helicity (or chirality). The equation we 
derived for spin states is 

{dt + V,, ■ V, + Fo.a,j5/« +sS.= -^ JCUk) - ICir\k) , (4.14) 

Jq TT L J relax 

where Si denotes all sources, both those from the fiow term p.5.135^|l.5.136jl and from the collision 
term ()2.58ti2759|l . while on the right hand side only the relaxation part of the collision term is left. It 
is now not difficult to break this equation up into four equations with s = ±1 and kz > and k^ < 0, 
respectively, and then combine those pieces with opposite s and kz to the equation 

{dt + vo^-V, + FoA.)Sh^ + hsign[k,]S,= / ^ /Co^,(A;)-/C^r"(A;) , (4.15) 

Jq Vr L J relax 

where Sfhi is the density of particles with helicity h minus the density of antiparticles with opposite 
helicity, precisely the quantity needed for the left-handed density ()4.13|) . Formally, this equation for 
helicity states has almost identical appearance as the one used in the heuristic WKB approach to 
baryogenesis (the classical force in the flow of the CP-violating density is usually omitted in these 
works), where the CP-violating source Si is approximated by the semiclassical force. What this 
analysis really shows is that the WKB approach, combined with a good guess work, which involves 
a correct calculation of the dispersion relation (obtained by transforming from canonical to kinetic 
momentum), and the right modus of insertion of the dispersion relation into the kinetic equation, 
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can lead to the correct CP-violating force. However, we emphasize that no work using the WKB 
approach has so far gotten the full source completely right. Another point is that the interpretation of 
h as helicity really makes sense only in the highly relativistic limit, or equivalently, in the symmetric 
phase, in which the tree level masses vanish. 

Furthermore, it is important to note that the explicit form of the source found in those works differs 
from our result: first, the contribution from the flow term to the source there contains only the 
semiclassical force 5Fi from ()1.5.137|) . but misses the other contributions we found. In section ^3 we 
make a numerical estimate of the importance of these terms. Second, the semiclassical WKB method 
is not in principle suitable to calculate the source that originates from the collision term. 

4.3. Fluid equations 

We now have the Boltzmann equations for the CP-violating distribution functions, and we know 
how to interpret them in terms of densities and currents and what is their role in the creation of 
a baryon asymmetry. A direct numerical solution of the Boltzmann equations ()1.9|1 fails due to 
their complexity, however, at least at present. So we have to rely on approximations to the full 
Boltzmann equations. A very convenient and with respect to the application to EWB quite popular 



approximation is a fluid Ansatz for the distribution functions 16 



for the fermions, and an analogous expression without spin indices for the scalar distributions. This 
form mimics the equilibrium distribution, but the chemical potential /i and the plasma velocity u, 
which are only functions of the spatial coordinates, allow for local fluctuations in the density and 
velocity of the plasma. In the given form the velocity perturbation accounts only for a net motion 
of particles in ^-direction. Because of the planar symmetry of the wall, this should be sufficient. 
The term including the wall velocity is due to the movement of the wall frame in the plasma rest 
frame. With the Ansatz ()4.16p the momentum dependence is explicit, such that, upon integrating 
the Boltzmann equations over momentum, only spatial dependencies remain. 

The CP-violating chemical potential and the velocity perturbation are caused by the interaction with 
the wall and hence are implicitly of first order in gradients. When we expand the fiuid Ansatz ()4.1fj|l . 
keeping only terms linear in /i, u and v^, and compare it with our decomposition p.5.128|) of the 
fermionic distribution function, leaving aside the CP-even correction, we can identify 

Sfsi± = -/?/oi(l - foi) [1 - Pvy,k^{l - 2/oi)] (-/isi± + k^Usi±) , (4.17) 

where /oi = l/(e^"'0' + 1), tuoi = (^^ + \md\fy/^. The form (fTTTjl of the CP-violating distribution 6f^ 
then suggests to define 

Hsi = Hsi+ - f^-si- (4.18) 

Usi = Usi+ -u_si-- (4.19) 

We insert the fiuid Ansatz into the Boltzmann equation for the fermionic CP-violating density 5/° 
and keep all terms up to second order in gradients and up to first order in Vyj. Then we take the 
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zeroth and first moment with respect to k^, that is we multiply by 1 and kz/ujQ, respectively, and 
integrate over the momenta. Since the flow term source is even in kz, it contributes to the zeroth 
moment equation: 



= - E {^^ts'^'^^s'^' + V^T'il,,Us'^') - J^ (/3r°>, + V^V'.^U,) . (4.20) 

s'i' i' 

In the following we do not consider the contribution p.5.136|) from the CP-even distribution function 
^/f™"^ to iSf°™, which in principle can be of the same order as p.5.135jl . The source from the collision 
term is odd in kz and therefore appears in the flrst moment equation: 

J {2nyuJo 

= Yl {^^(3^lU'^'s'^' + rL>v«s',') - J2 (^-/^rj>, + ri>,) . (4.21) 

s'i' i' 

Here dots denote the derivative with respect to time, primes are 2;-derivatives, and we introduced the 
dimensionless integrals 



W = P'^'''' l^..{l-h.)^^ (4.22) 



Jai,,i = /3'+'^-' / 7Tr^/o.(l - /o.)(l - 2/0.)^ . (4.23) 



With the fluid Unsafe the momentum dependence of the 6fg and 6^^ is specifled, so that the integrals 
in the relaxation term (j3.17|) can be performed, leading to the relaxation rates r^j^s'i' and T^j, which 
connect different fermionic species with each other and fermionic particles with scalar particles, 
respectively. Some explicit expressions for the rates calculated in our one-loop approximation can be 
found in appendix |D1 Note however that with our conventions the above fluid equations imply that 
the actual rates, by which the chemical potential and the fluid velocity perturbations are attenuated, 
are given by ^ Ji sn' / ^oo,i and T^^ ^^i / l2i^i., respectively, which are about one order of magnitude larger 
than Tll^,^, and T]lg,^,. 

The equations for the scalar particles can be obtained from the ones for fermions in a simple way: 
one has to replace \m\^ by M^, remove the source from the flow term, replace the coUisional source by 
the scalar one, use the appropriate relaxation rates for scalars, and replace the coefficients / and J by 
their scalar counterparts, whereby the Bose-Einstein distribution replaces the Fermi-Dirac function. 

The fluid equations ()4.2Uffl?^T|) are well suited for a numerical treatment. They consist of a system 
of flrst order differential equations, which is a comparatively simple problem. Note that only the 
time- and the z-derivative have survived. This is a consequence of the fluid Ansatz, which in the 
above form effectively allows for spatial variations only in the ^-direction. Since in the wall frame 
the mass depends only on the z-coordinate, the problem is stationary in essence and we could drop 
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the time derivative. We argued in the beginning of section 11.5.1.21 that keeping the time derivative 
could allow a treatment of non-equilibrium initial conditions. But in order not to violate the spin, 
which is conserved by the symmetry of the problem, these initial conditions have to have quite a 
special form, given by equation p.5.23|) . so that their physical relevance is at best limited. Keeping 
the time derivative is however a convenient tool to solve the equations. At a first sight it seems to 
be simpler to have only ordinary differential equations instead of partial ones. But to solve these 
ordinary equations is actually quite tricky, since boundary conditions at infinity have to be satisfied: 
at sufficiently large distances apart from the wall we expect the system to be in chemical equilibrium. 
Starting with some initial values at one side of the wall usually does not lead to the correct chemical 
equilibrium on the other side, so that the initial value has to be fine-tuned. A simpler procedure 
is to keep the time derivative and thus simulate relaxation of the system toward a stationary state, 
starting from some well-chosen initial functions fio{z) and Uq{z). Because of the interplay between 
the dissipative effects in the relaxation term on the one hand and the sources on the other, this 
procedure leads to stable stationary solutions relatively fast. 

Once the solutions of the fiuid equations are known, the axial density can be calculated. Insert- 
ing ()4.17p into ()4.1H) leads to the relation 

f = -J2'J-(^f^fo(^-fo)f^[^s + v^P{l-2fo)f,s). (4.24) 

It is interesting that the main contribution to the axial density here comes from the plasma velocity, 
since the chemical potential is suppressed by an extra factor v^. 

4.4. Solving the fluid equations 

We now solve the fiuid equations in a toy model consisting of one fermionic particle species. Despite 
of the marked simplicity of this model, the explicit solution of the equations will be instructive, 
and in particular will help to answer three questions. First, we would like to know how big is the 
collisional source ()2.59p . when compared to the fiow term source ()1.5.135|) . Since these sources appear 
in different equations, a direct comparison is not possible. But we can study the influence of the 
presence or absence of the collisional source on the axial density. This should allow us to judge if this 
contribution could be important in an actual calculation of EWB in a realistic model. Second, we 
found that previous work based on the WKB method missed a part of the flow term source p.5.135|) . 
Are the new terms relevant? And flnally, we have seen that in the highly relativistic limit on the 
level of Boltzmann equations the use of spin states or helicity states is equivalent. As soon as an 
approximation to these equations is made, however, as for example the fluid Ansatz, differences occur. 
We can solve the fluid equations in both pictures, compute the axial density and compare the results, 
in order to study by how much the equivalence between spin and helicity states is broken on the level 
of fluid equations. 

The bubble wall can be modeled by a hyperbolic tangent 

|m| =mo-(l -tanh-^) , ^ = ^077 (l - tanh-^) (4.25) 

for both the absolute value and the phase of the fermionic mass. Here L^ characterizes the wall 
thickness. In contrast to the fermionic mass parameter mo, which influences the form of the source 
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FIG. 9: The axial density generated by the flow term source, the "old" flow term source, and the collisional 
source, respectively, for different values of the fermionic mass parameter. The pictures on the left hand side, 
computed with the rates 1)4. 26() . show no diffusion, while in the pictures on the right, computed with 1)4. 27() . 
the effects of diffusion can be seen. For this choice of parameters the contribution of the collisional source 
is negligible. The difference between the results for the two flow term sources is apparent. The bubble wall 
with width L^ = 10/T sits at z = 0, the broken phase is at the left. Note the different scales of the z-axes. 



34 





FIG. 10: The effects of diffusion can be nicely seen here, where we plotted the axial density produced by 
the same source, but with different sets of relaxation rates, (|4.2Hj) and (|4.27j) . For comparison we also show 
the flow source. The axial density computed with the larger relaxation rates is not transported into the 
symmetric phase, but sits right on top of the source, while for the small rates a diffusion tail is clearly visible. 
For small masses diffusion is more efficient. 

and therefore also the form of the resuh, 6q is simply a common pref actor to all source terms and 
therefore only scales the result of the computation. We choose 6q = 0.1. The width of the bubble wall 
is set to Lyj = 10/T, a typical value for the MSSM, and the wall velocity is chosen to be Vw = 0.1. 
We want to investigate the collisional source for the fermionic equations, so the system must also 
contain scalar particles. But since we are interested only in the dynamics of the fermionic particles, 
we assume a thermal distribution for the scalars. Both fermions and scalars get their masses from 
the interaction with the Higgs field, and therefore the scalar mass is proportional to \m\. To keep 
things simple, we set all relaxation rates in ()4.20|) and ()4.21|) to zero, except of F^^ and Fj". For 
these rates we use two different sets of values, 



FL^ = T/25, Tl: = T/20, s = ±l 



(4.26) 



which are so big that there will essentially be no diffusion of particles in front of the wall, and 

F°^ = r/120 , F,^^ = T/lOO , s = ±1 , (4.27) 

which are small enough to have a notable diffusion. Note that both fluid equations contain the time 
derivative of the chemical potential and of the plasma velocity. We flrst build linear combinations of 
these two equations, such that we obtain one equation which contains only the time derivative of /x, 
and another equation in which only the time derivative of u appears. The temporal evolution of these 
equations can then be simulated numerically in a straightforward manner, where we choose fi and u 
to be identical to zero as initial values. We emphasize again that this evolution is a physical temporal 
evolution only for the very special initial conditions p.5.23j) . But in any case it is a convenient way of 
obtaining the solutions of the ordinary differential equations in z, which emerge by setting the time 
derivatives to zero, without having to worry about boundary conditions. 
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FIG. 11: The axial density computed by solving the fluid equations for spin states and helicity states, 
respectively. The results are very similar, but for both values of the fermion mass there is a clearly visible 
difference. 

To address the first two questions we mentioned above, we solve the fluid equations for three different 
cases: with the flow term source, with the "old" flow term source, that means using only the first 
term of p.5.135j) . as it was done in the WKB calculations, and with the collisional source only. 
The computation is done for three different values of the fermionic mass parameter mo, and for the 
collisional source we set M = 3.0|?7i| and use |?/p = 0.1 for the coupling constant. In figure IHl the 
resulting axial density is shown as a function of the z-coordinate. On the left hand side are the 
results obtained by using the larger relaxation rates ()4.2fi|l . while for the pictures on the right hand 
side the smaller rates ()4.27|1 have been used. For the parameters chosen, the collisional source has 
no notable effect. Repeating the calculation with different wall velocities and scalar masses does not 
increase the collisional contribution significantly. Note that for smaller relaxation rates not only a 
diffusion tail develops, but the amplitude of the axial density grows, too. 

The axial density originating in the "old" flow term source overestimates the one from the true flow 
term source, with a difference increasing with the mass parameter. While for small mass parameters 
mo < T both sources give about the same result, the difference becomes significant for larger masses. 
We conclude that a computation of EWB with the correct flow term source will most probably lead 
to different results for the baryon density than the WKB calculations. 

In order to illustrate the difference between the situations with and without effective diffusion, we 
show in figure HHl the axial densities calculated with the two sets of rates ()4.2(i|l and ()4.27j) in direct 
comparison. 

Working with spin states instead of helicity states leads to two important changes on the level of the 
fluid Ansatz. Because in the Boltzmann equation ()4.15|) the source is multiplied by sign{kz), the flow 
term source for helicity states will appear in the flrst moment equation, and the collisional source 
will go into the zeroth moment equation, just opposite to the fluid equations for spin states. The 
other difference concerns the axial density, which in terms of chemical potential and plasma velocity 
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now reads 



f 



J2^J T^^^o^^ - ^o) [f^h + t'^./Sll - 2fo)u^,) . (4.28) 

Compared to ()4.24p for spin states, the roles of chemical potential and plasma velocity have changed, 
in a way. The question is now, whether the effects of both of these changes cancel each other. So we 
solve the fluid equations also for helicity states and show the result for the axial density in figure ITTl 
in direct comparison with the density obtained by working with spin states. One can see that the 
results differ by about ten percent. For small masses the fiuid equations for spin states produce a 
larger axial density, for high masses the equations for helicity states produce the larger result. Since 
both cases are an approximation to the full Boltzmann equations, this difference can to some extent 
also be regarded as an estimate of the error one makes by going from Boltzmann to fluid equations. 

4.5. Diffusion equation and sources 

Here we make a rough analytic comparison between the flow and collision term sources. In the 
stationary limit and neglecting the terms of order f ^ <^ 1 (as usually, in the wall frame), Eqs. ()4.2()l - 
I4.21|) can be recast as 

I2lX^ - /om/5^^«. - t;. JoM/3=^^/is. -(3'f 7^«5f°- (4.29) 

2 2 J [Zn}-^ 



s'r 



E [rS.v/?/^s'.' --.(riv-p^r°i:,,,)«.,,] +Y. \^Z'f^f^^'--^(^l^'-j^^^:i' 



h2,^(3|^s^ + 'Vw[^7J ^22,^/5 \md\i Usi - (3 / — — ^ sS^ (4.30) 

s'i' ' i' ' 

where we kept only those sources which are not suppressed by the wall velocity. (Recall that the 
sources iSf°™ and 5?°'' are both linear in f^.) From Eq. (J4.3UJ) we can easily write an expression for 
the fluid velocity. 



1 



-plu 1 „, fj Ioi,iJ42,i \ Q2\m ,|2^ 
^ si,si + ^w\J22,i 21^^ )P \^d\i 






-4-7, '^'^^•* R^ f _^!^. Cflow _ /Q3 / ^^^ ^z <.coll 
+"-7^^ J (2^0^^'^^ ^ J J2^^W^'^^ 






,(4.31) 



and similarly for /Xsj. Upon inserting Eq. ()4.31|) and its derivative into the continuity equation ()4.29|) . 
we get a rather complicated second order diffusion equation for the chemical potential, from which 
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FIG. 12: The rate F^" as a function of the fermion mass |rTT-|. The scalar mass M is a multiple of |?tt.|. At 
the ratio M/\m\ increases, the rate peaks at an asymptotic value Fj" — > (0.9 — 1.0) x lO^^T. The spin-flip 
rate F^'^^ is negative and smaller by a factor four, approximately. 



we list just a few relevant terms, 

'22.i J-22.iJ-21'i 



n // 7~i \^ f si,s'i' Ji2,i j^Oii \n I 



(4.32) 



D 



SI n2 
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(jb rv fXi y 



■ss: 



coll' 



T I \2' 



l22,i' J (27r)3c^or ' ho,i' 2 

where we defined a diffusion coefficient 



l2^^ 



■ssr'' - 



'00,j 
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(2^ 
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0, 



hl,ih2,i 



'00,4 
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si, si 



+ 0{v^) 



(4.33) 



We can neglect the contribution of the coUisional source term oc \md\f , unless the diffusion constant 
is very small and Irridli ^ T. From Eq. ()4.33p it is then clear that as the diffusion parameter Dsi 
grows (or equivalently, when Tl^g^ becomes small), the collisional source becomes more and more 
important. To give a rough quantitative comparison between the sources, we first note that the flow 
and the collision term sources p.5.138|l . ()2.()()j) can be approximated by 



P' 



P' 



(2^'^^ ~ 2x10 — (^-J, + — T, 



{2lT)^Uoi 



(4.34) 
(4.35) 



where we made use of dz ~ l/L = T/10, 6 r^ Oq = 0.1, \md\i ~ mo, and Vw = 0.1. Now from figures lL5l 
and [7| we see that X^ and X^ peak at about one, and are typically one order of magnitude larger than 
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T, and we conclude that the flow term source (J4.34J) is typically about two orders of magnitude larger 
than the collision term source ()4.H5|1 . When viewed as sources in the diffusion equation ()4.H2|1 . we 
first observe that both sources appear at second order in gradients, and are hence suppressed as 
1/L^, so changing the wall thickness L^ should not affect much the relative strength of the sources. 
Analogously, the relative strength of the sources is not changed by changing the wall velocity f^, 
Yukawa coupling \y\, or CP-violating phase 9. Second, making a comparison of the collisional and 
flow term sources in Eq. (j4.32p at a maximum |md|j ~ T leads to the conclusion that the sources are 
comparable when the diffusion constant becomes Dgi ~ 300T~^, such that T\^^^ ~ 0.03/^^^ ~ lO^'^T 
is the rate below which the collisional source should dominate j^. A self consistent quantitative 
comparison of the two sources can be inferred from the above analysis and Figure ^1 where we show 
the rate T],^^^ as a function of the masses \m(i\i and Mj. In the parameter range where the sources 
peak, \mci\i — T and Mj ~ (3 — 4)T, the rate Fj^^j ~ (3 — 5) x lO^'^T. This then implies that, when 
compared with the flow term source, in the region of parameter space where both sources peak, the 
collisional source is a few times smaller. This collisional source suppression can be attributed to a 
phase space suppression of 5,^°'', to which only particles with relatively small momenta contribute, 
which are not so numerous. We expect that a similar conclusion can be reached when higher order 
loop processes are considered, both as regards to the collisional source as well as diffusion constant 
contributions. 
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5. CONCLUSION 

With the present paper we complete our work on a controlled first principle derivation of transport 
equations for chiral fermions Yukawa-coupled to complex scalars by a discussion of the collision term, 
whereby we approximate the self-energies by the one-loop expressions. 

We calculate, from first principles, the CP-violating source appearing in the collision term of the 
fermionic kinetic equation. This represents a first calculation in which both semiclassical and spon- 
taneous sources are obtained within one formalism. We also present a first computation of a CP- 
violating source from the scalar collision term, which arises from nonlocal fermion loop Wigner 
functions, and can be traced to the CP-violating split in the dispersion relation for fermions. Next, 
we derive the fiuid equations associated to the Boltzmann equations found in Paper I and solve them 
(by a novel method) for a simple case, including only fermionic particles, which allows us to make a 
comparison between different sources. It turns out that the collisional source is subdominant when 
compared with the semiclassical source coming from the semiclassical force in the flow term of the 
fermionic kinetic equation, unless the fermions diffuse very efficiently. 

We discuss the relation between the helicity and the spin quasiparticle states and show that in a 
highly relativistic limit the Boltzmann equation for spin quasiparticles reduces to the one for helicity 
quasiparticles. Fluid equations are, however, not the same. In particular, we show that in the 
fiuid equations for spin quasiparticles the semiclassical force sources the continuity equation, while 
the collision term sources the Euler equation; exactly the opposite is true for the fiuid equations of 
helicity quasiparticles. 

The question of thermalization rates in the collision term of the Kadanoff-Baym equations is also 
addressed. We calculate the rates with the one-loop self-energies, for which only scalar absorption 
and emission processes are taken into account. We point out at the shortcomings of these rates: 
as a consequence of the energy momentum conservation enforced in the on-shell limit, the rates are 
suppressed for small masses, and vanish in the massless limit. By establishing a relation to more 
standard rates calculated with helicity states, we outline how one could approximate scattering rate 
calculations from the Kadanoff-Baym equations with two-loop self-energies, which are needed for the 
derivation of quantitatively correct fluid equations. 

We conclude by making a few comments on further directions of the research presented here and 
in Paper I. Solving the Kadanoff-Baym equations with the self-energies truncated at two-loop or- 
der, needed for modeling baryogenesis at a flrst order electroweak phase transition, would require a 
signiflcantly more effort than what was needed in this work, and thus would be justifled only when 
quantitative consistency checks of a particle physics model would be needed. Testing the accuracy of 
the fluid equations based on the Kadanoff-Baym equations with one-loop self-energies is feasible, and 
should be performed. Furthermore, it would be of interest to develop a basis independent formalism 
for mixing fermions valid to order /i in a derivative expansion, and in particular to study the signif- 
icance of the mixing in the situations when mass eigenvalues are (nearly) degenerate. Moreover, an 
inclusion of the hermitean self-energies would be also desirable. Apart from providing thermal mass 
corrections, these self-energies mix spin and may provide an additional source of CP violation (which 
we expect to be a of a similar strength as the collisional source). Finally, when fully developed, this 
formalism should be used for quantitative baryogenesis calculations. Important examples comprise 
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baryogenesis mediated by charginos or neutralinos of super symmetric extensions of the Minimal Stan- 
dard Model, baryogenesis mediated by (top) quarks or scalars in two Higgs doublet models, etc. The 
interest in these type of calculations is particularly well motivated in the light of the next generation 
of accelerator experiments (LHC, NLC), which will provide a new experimental input on the various 
parameters of the electroweak scale physics. Moreover, novel experiments on electric dipole moments 
will probe in more depth CP-violation in the electroweak sector. 
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APPENDIX A: CP- VIOLATING SOURCES IN THE SCALAR COLLISION TERM 

In this Appendix we show how to perform analytic evaluations of {seven out of eight integrals in) the 
integral ^TMi 

Recall that this integral is important because it yields a CP-violating source in the scalar collision 

termC^o ^HIB- 

We argued in section IZJTl that, when evaluating (jAljl . it suffices to work with the thermal equilibrium 

particle densities ()2.14j) . which we rewrite as 

^feq = ^- [^(^0 - ^0)/0 + 5{h + C^0)(1 - /C 



^feq = -^- ['^(^o - cuo)(l - /o) + S{k, + cuo)/o] (A2) 



where ujq = y k"^ + \m\'^ and we used the following equilibrium definitions for the particle and an- 
tiparticle distribution functions 

/o = '^eq('^o) = —^ (plasma frame particle distribution function) 



1 



/o = 1 — neq{—uJo) = ^^^^^ ^ ^ (plasma frame antiparticle distribution function) . (A3) 



Note that in thermal equilibrium in plasma frame and to leading order in gradients we have, /o = /o- 
We can now easily evaluate the fcg ^^d feg integrals in XTJ^. The result is 



I,.ik) . X> (.) ^ -1 / '-^Srk + k' - k")^^ 

S{ko + UJ',- .;^')/o(l - /o ) + '^(^0 - ^0 + ^o)(l - /o)/o' 



lo'qUj'^ - k[ ■ k'' 






0^0 

u'u'' + k' ■ k'! 



+ 






5(A;o + 4 + ^o)/o/o' + '5(^o-^;-^^')(l-/o)(l-/o) }, (A4) 



where u;^ = ^J k' + \m\^, u'f^ = \J k" + |m|2, i;^ = \jk'^ \ jmp, and /q = rao(c^o)' /o = ^^0(^0)5 
/g = 1 — 'n,o(— cUq), /q' = 1— no(— cJq). The former two terms in ()A4|) are contributions from absorption 
and emission of a scalar particle. The latter two arise from annihilation of a scalar particle with a 
fermion-antifermion pair and the fermion-antifermion pair creation from a scalar particle, respectively. 
These four processes are illustrated in figure |H1 
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Upon performing the /^''-integral, Eq. (jA4j) can be reduced to the following integral 



1<t>o{k) 



1 



'3l' 



d'k 



X 



167r2 ./ {Cj'^y 

K^ + h-K 



^0^0 



+ 1 + 






5{k, + J,- ^^')/o(l - /o ) + 5(^0 - ^0 + ^o)(l - /o)/o' 



<5(fco + a;[, + a;[,')/o/o' + 5(^0-^^- ^(,0(1 -/o)(l-/o) , (A5) 



where now uj'q = a/Zc'^ + |mP, cUq = a/A;^ + 
convenient to break the integral 



19 '/ 



(fc + A?')^ + I^P, ^[i ■ Ml = K\' + ^11 ■ Ml- It 



IS 



oo /-oo />27r 

3;^' / j;,/ / iJ jiJ 



d'k 



oo Jo 



dk[ / fcurfA; 



(A6) 



where k\\ ■ k'u = k\\k',, cos (p' . We can perform the angular 0-integral in ()A5jl with the help of the 
(5-functions. Note first that 

6{±ko + cu',-u'^) = 9iTko)9{k'k''-{Po-kXy) 



X 



Xn = COS I 



'^o — ^0 ~ ^<t>^ 



II II V "^0 

(MV2) - u^u'o - fc.fc^ 



5(0'-0'o) + 5(0' + 0[,) 



hk' 



(AT) 



such that in the first (5-function only the antiparticle pole, ^(— /cq), contributes, while in the second 
(5-function only the particle pole, 9{ko), contributes. Similarly, we have 



6{±ko + u;', + u;'^) = e{Tko)e{k''k'^ - {(3^ - k^f) 



X 



k\\k'>/T^ 



X' 



l2 



0[,)+5(0' + 0'o) 



i^O ~ ^(/> ^ ^0 ' 



(A8) 



The product of ^-functions in (IA7IIA8|) stems from the requirement that |xg| < 1, and can be conve- 
niently expressed as 



M > 2\m\ 
kk' > \(3o\ , (^0 = ^ ^</.^o • 



K, <K< M„ 



(A9) 
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The result of integration over the (5-functions is then 
2^00 (fc) = 






e{u^ - u'o) [e{-ko)fUi; + e{ko){i - /^)(i - fd)] 



^0 —'^tl>~'^0 



(AlO) 



To evaluate the fc^'i^^^grals in ()A10|) the following integrals are useful. They can for example be 
performed by deforming the contour of complex integration - the closed curve that contains the 
poles k'^ = ^22 1 ^ i^to ^ contour around the complex poles k'^ = ±i|m|, and evaluating the sum of 
the associated residues: 



Xi 
Xs 



dk' 



■K 



^- 



1 



K, K + \M'^{K~K){K-K) H ^{\m\'-KK) + AMiK + K) 



k^dk^ 



1 



dk' 



-S 



TT 



ViH''K,k'J + ^rn\iK, + KJ 



where 



KdK 

k'^kl + {131 - k^k'^) 



TT 



Tik, 



PO: 



k' k' 



kr, + k^. 



2kJo 



^2 ' zi • 22 ^2 



The integrals Iq and Xi are generally expressed in terms of complex functions of the form 

1 k 1 



(All) 



(A12) 



(A13) 



^{\m\^-ki^k'J+t\m\{k',^ + k'J k'^{k,-c,){c^-k,) 
where c_ and c+ are the (complex) roots of the equation (|mp — k'^^k'^^) +i\m\{k'^^ + /c^J = 0, that is 



c± = 

We can thus reduce ()A10|) to the form 

1 



i\m\(3o uJq 



k'^ ±p^/^^^''-/5. 



(A14) 
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_2\m\ ^{k,-c.){c+-k,) ^{k,-cJ){c+-K) 

e{Jo - ^^) [0{-ko)m - /o ) + 0{ko){l - fo)fo] 

- e{u^ - u'o) [^(-fco)/o/o + ^(^o)(l - /o)(l - /o')] 



k 



0— '^o~ 



U),-, —Ul,,—LxJ^ 



'o=^4--^o 



(A15) 



The final fc'-integral cannot be performed analytically. In section 1^.1. II we use this integral expression 
to compute the CP-violating collisional source in the scalar kinetic equation. 
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APPENDIX B: CONSISTENCY OF THE FERMIONIC KINETIC EQUATIONS 

In section 1X31 we have presented a derivation of the fermionic kinetic equations. As a consequence of 
spin conservation, the Chfford algebra of 16 functions is reduced to four (a priori different) functions 
gl^ (for each spin s = ±1). Furthermore, the constraint equations p.5.56MT.5.59|) relate different g^, 
resulting in one independent function, for which the standard choice is the vector particle density on 
phase space, g^. Here we prove that all kinetic equations p.5.80III.5.83|) are equivalent both at the 
level of flow, as well as collision, terms, implying consistency of our kinetic equations. First we show 
the equivalence of the flow terms and then we consider the collision terms. 

We shall not give a proof for all three equations, but instead restrict ourselves to the kinetic equa- 
tion ()1.5.83|) for gl: 

Uhdt + ^11 ■ d)g', + sd.g^, + 2(m, - \m:dl)g{ - 2(m, - \mldlX2 = ^3 • (Bl) 

We insert the expressions ()I.5.62|) - p.5.63|) for gl and (7I (accurate to second order in gradients, since 
the kinetic equation for g^ contains zeroth order terms), to get 

i (sK + ^M^e'du}] ^{kodt + k\\ ■ d)g', (B2) 

ko ko 

where K.^, C^ and C^ are defined in Eqs. p.5.84|) and p.5.60|) . respectively. Now we take the derivative 
of the constraint equation ()I.5.68|) with respect to z (the coUisional contribution to this term can 
then be neglected), and we find 

-\mfg^, + ^{H'ejg^, + (k' - \m\' + f |m| v) g^, = . (B3) 

ko V ko J 

Next, we use this relation to replace the term —\m\'^gQin the second line of ()B2|) . After some algebra 
this leads to 



ko\ 2ko y[ko 2ko 2kr 



,^ = /C^-^q + ^C^ (B4) 

We recognize this equation as some operator times the left hand side of the familiar kinetic equation 
for g^, thereby proving that the flow term of the kinetic equation for g^ is consistent with the one for 
gQ, which thus carries no new physical information. The analysis of Eqs. p. 5. 81^1 and p.5.82|l for gl 
and (^1 is somewhat more complicated, but the final conclusion is the same. 

From Eqs. p.5.80|) and ()B4|) it follows that one of the conditions required for the consistency of the 
collisional sources is 

^3 - l^W + -l-<-2 = -7-1^0 ■ (B5) 

fco Ko Ko 
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In oder to show that this relation is indeed satisfied, we make use of the zeroth order coUisional 
source ()2.54j) and the first order coUisional source ()2.5(ij) . whose precise form is 



7;M /■ rl^h'rl^h" ^^ 



4 '" "^^"'y (27r)4 
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We find 



/Cq 


= ^^ 
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= f^k 


Cl 


= ^l 
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= ^l 






\/tQ ftp 

ss \m\ U ^7—^3 h rrii^^^- — h ss rnf^^^ 

Kq Kq Kq Kq ko_ 



ko k'^ 



where /C|[] and /C| [] denote the functional 



/I \'2nl^^h ^ . I kz k. . . ^ 

SS \m\ U ^7—^3 h m^^^- — h ss m^^^ 
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(B7) 



}Cl[-]^^(3^^v. 



K [•] = '-^hv,vr 



l/r .. .. fd^k'd^k" ,4,, , , , , „, ^ [^ , ,A;oA;^ - A;,, ■ k[ ^ ^^ 



(2vr)4 
d^k'd^k" 



(54(A;-A;' + A;")$^(l + ss' 
5\k-k' + k")Y,{^ + ss' 



kok'o 

f '-' U II 

kok'o 



^oe<(fc)^o4(A:')^A<(A:")[-] 

fB8l 



From this it immediately follows that the consistency condition ()B5|1 is satisfied. Note that the last 
two terms in each of the sources Cf and C| originate from the first order coUisional source ()B6|) . 

By considering the kinetic equations for gl and (yfl, one can show that similar consistency relations 
hold, which completes the consistency proof for the fermionic kinetic equations. 
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APPENDIX C: CP- VIOLATING SOURCES IN THE FERMIONIC COLLISION TERM 

In this appendix we show the details of the calculation of the fermionic collisional source term, which 
we omitted in the main text. We show here the calculation in the mixing case, the non-mixing case 
is simply obtained by setting \m\f = |?t2|^ = |mp. We begin with equation ()2.70|) : 



j (2^TZ7'^^''°' ^ ''°^' ~ ''^)^o^4(i + /o 



s[ \m 






The k' integration is performed in cylindrical coordinates, d^k' = dk'M,dk'nd(f), where is the angle 
between /cy and ki, , and k\\ and k',, are the absolute values of the parallel momenta, respectively. We 
substitute 

x = k-k' = kiik\^ cos (I) (C2) 

and instead of integrating over the absolute value of the parallel momentum we integrate over k' = \k'\: 

d^k' = / dk',- / dk'^e{k'^ - kl) / , d{h\k[ - \x\) . (C3) 

^-oo 2 Jo y_oo /(A;||A;[|)2-x2 

Since u'^ = {{k — k'Y + M^) ^ = (A;^ + A;' — 2x — 2kzk'^ + M^) ^ is a function of x, we can rewrite 
the remaining (5-function 

5{u:o + u:'o- u'l) = ul5{x- (/?o - k^)) , (C4) 

where we used the abbreviation 

/3o = ^ ^ ^oi^oi • (C5) 

With the 5-function the x-integration can now be performed trivially. With a not so trivial calculation 
one can show that for arbitrary /5o the relation 

e{k\\k\ - (/?o - k,k';))6[k'' - k'^) = e[ek" - ^l)e[k'^^ - k',)e{k', - k'^J (ce) 

holds, where k'^^^ are the roots of (A;||A;[j)2 - (/3o - k.k'^f = 0, see (jMI)- The function 9{k'^k''^ - (3l) 
contains a mass-threshold. It can be found by trying to determine those values of k' for which the 
6'-function is non-vanishing: 

Pi < Pk'^ ^ -kk' <(3o< kk' . (C7) 

The latter two conditions can only be satisfied if the quadratic equations 

fc'Vil' T k'k{M^ - \mi\^ - |m,f )2 (C8) 

i;2i l2ii |2| |2 ^ / ^4-2 I |2 I |2\2 n 

+ k \mj\ + \mi\ \mj\ [M — \mi\ — \mj\ ) = , 
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have real solutions, which is only the case if 



M^ > |mip+ |m,-p 



(C9) 



Now we can write the source as 
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For the fluid equations the integral over all momenta is needed. With the observation that the 
integrand is odd under k^ -^ —k^ and k[^ -^ ~K^ ^^ can immediately write 



d^k 



(27r 



/Cn ,■ = . 



»4 0:« 



(cii: 



where the index ± denotes the integral over positive and negative frequencies, respectively. The 
integral with an additional factor k^/uQi can be evaluated further by making use of the integrals ()A11|) : 
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(C12) 
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APPENDIX D: COLLISIONAL RATES 



With the fluid Ansatz (j4.16|) the momentum dependence of 5g is completely specified. 

The evaluation of the integrals runs along the same lines as for the coUisional source in the first part 

of this appendix. We just give the results. The contribution to the zeroth moment equation is 
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In the first moment equation the collision term is 
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where, besides the source term, there are the rates 
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